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PREFACE 
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(AFSC),  under  Program  Element  65807F.  The  results  of  the  research 
were  obtained  by  ARO,  Inc.  (a  subsidiary  of  Sverdrup  &Parcel  and 
Associates,  Inc.),  contract  operator  of  A EDC,  AFSC,  Arnold  Air 
Force  Station,  Tennessee,  under  ARO  Project  Numbers  RF217  and 
RF414.  The  author  of  this  report  was  John  M.  Pelton,  ARO,  Inc. 

Data  analysis  was  completed  on  June  30,  1974,  and  the  manuscript 
(ARO  Control  No.  ARO- ETF-TR-75-92)  was  submitted  for  publica- 
tion on  June  25,  1975. 
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1.0  INTRODUCTION 


At  the  present  time,  the  gas  turbine  engine  industry  has  no 
generally  accepted  set  of  thermodynamic  properties  for  the  constitu- 
ents of  the  working  fluids  or  method  for  making  thermodynamic  calcu- 
lations. Therefore,  each  segment  of  the  industry  from  manufacturer 
through  user  may  use  slightly  different  basic  thermodynamic  data  and 
calculation  methods  to  arrive  at  the  thermodynamic  properties.  The 
present  condition  is  similar  to  that  which  existed  in  the  rocket  motor 
industry  prior  to  the  work  of  the  ICRPG  (Ref.  1)  in  establishing 
methods  for  performance  calculations  and  data  reduction. 

As  the  operating  environment  of  gas  turbine  engines  has  gone  .to 
higher  temperatures  and  Mach  Numbers,  the  differences  in  perform- 
ance of  competitive  engines  are  often  quite  small.  In  fact,  some  per- 
formance differences  are  of  the  same  order  of  magnitude  as  the 
change  that  may  be  expected  from  using  different  sets  of  thermo- 
dynamic properties  of  the  constituents  in  the  working  fluids.  Also 
because  of  the  increased  interaction  between  the  engine  and  other 
components  of  the  airframe,  the  need  for  a consistent  method  of 
arriving  at  thermodynamic  properties  becomes  even  more  important. 
The  increased  use  of  mathematical  models  for  describing  engine  per- 
formance has  also  increased  the  need  for  the  use  of  consistent  thermo- 
dynamic properties  of  the  working  fluids. 

The  gas  turbine  industry  has  developed  from  industries  with 
diverse  backgrounds  and  each  segment  has  developed  calculation  pro- 
cedures based  on  its  own  past  history  and  present  needs.  The  appear- 
ance of  the  computer  as  a routine  tool  for  turbine  engine  analysis  in 
even  the  smallest  operation  now  provides  the  possibility  for  use  of  a 
consistent  set  of  thermodynamic  properties  by  all  segments  of  industry 
and  government  if  the  properties  can  be  coded  for  use  with  small  com- 
puter memory  and  calculation  times. 

In  an  effort  to  develop  a consistent  thermodynamic  model,  certain 
guidelines  were  set  up  as  follows: 

1.  Accessible  constituents  and  constituent  properties  that 
are  readily  located  so  that  orderly  updating  is  possible 
as  compositions  for  standard  air  and  fuel  change  or 
thermodynamic  property  data  are  improved, 
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2.  Acceptable  memory  requirements  and  calculation 
times  for  use  in  on-line  data  reduction  systems,  and 

3.  Contains  real  gas  effects. 

Certain  items  are  required  in  such  a consistent  model,  and  these  are 
identified  below: 

1.  The  composition  of  air, 

2.  The  composition  of  hydrocarbon  type  fuels, 

3.  A set  of  thermodynamic  properties  for  the  constituents 
of  working  fluids,  including  real  gas  effects,  and 

4.  A method  of  computing  the  thermodynamic  properties 
of  air  and  exhaust  gas,  including  real  gas  effects. 

Each  of  the  above  items  was  examined  by  reviewing  the  available 
literature  and  the  present  turbine  industry  usage  in  conjunction  with 
the  available  computer  technology  for  on-line  data  reduction  calcu- 
lations. 

The  pressure-temperature  regime  for  air  as  a working  fluid  is 
shown  in  Fig.  1.  The  regime  is  from  0.  1 to  600  psia  and  from  300 
to  1,800°R.  These  limits  were  determined  from  current  and  projected 
future  requirements  for  air  as  a working  fluid  in  turbine  engines.  A 
number  of  air  compositions  were  reviewed.  They  include  the  work  of 
Self  (see  Edmunds,  Ref.  2)  and  Keenan  and  Kaye  (Ref.  3),  Hilsenrath, 
et  al.  (Ref.  4),  Touloukian  (Ref.  5),  Banes  (Ref.  6),  Brahinsky  and 
Neal  (Ref.  7),  Brown  and  Warlick  (Ref.  8),  and  the  U.S.  Standard 
Atmosphere  (Ref.  9). 

The  composition  of  the  hydrocarbon  fuel  was  defined  as  CnHxn 
where  "x"  would  be  determined  by  the  particular  grade  fuel  in  use. 

There  is  a very  limited  amount  of  thermodynamic  data  from  a 
single  source  that  covers  the  necessary  pressure  and  temperature 
regime  for  the  aiT*  and  exhaust  gas  constituents.  The  regime  for  air 
has  been  shown  in  Fig.  1,  and  the  regime  for  exhaust  gas  is  shown  in 
Fig.  2.  The  regime  for  the  exhaust  gas  is  from  0. 1 to  600  psia  and 
from  600  to  4,000°R.  The  fuel-to-air  ratio  (f)  to  be  considered  will 
be  from  zero  to  stoichiometric.  These  model  limits  were  based  on 
current  and  projected  future  requirements  for  calculating  exhaust 
gas  thermodynamic  properties  in  the  turbine  engine  industry.  Because 
of  the  overlapping  of  air  and  exhaust  gas  constituents  such  as  nitrogen. 
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oxygen,  carbon  dioxide,  argon,  and  neon,  the  temperature  range  for 
these  constituents  will  be  from  300  to  4,000°R.  Constituents  such  as 
water  vapor  require  thermodynamic  properties  from  600  to  4,000°R, 
while  the  products  of  chemical  dissociation  require  thermodynamic 
properties  from  approximately  2,000  to  4,000°R.  The  two  most  com- 
plete sources  of  properties  are  JANAF  (Ref.  10)  and  McBride,  et  al. 
(Ref.  11).  Other  sources  reviewed  included  Hilsenrath  (Ref.  4), 
Touloukian  (Ref.  5),  Din  (Ref.  12),  and  Stull  (Ref.  13).  None  of  these 
sources  was  complete  in  terms  of  required  constituents,  temperature 
range  or  recent  review  and  updating,  but  the  JANAF  data  (Ref.  10) 
were  chosen  because  of  the  temperature  range  and  the  on-going  work 
in  up-dating  the  data  as  more  information  becomes  available.  Two 
constituents,  argon  and  neon,  are  not  included  in  these  tables;  there- 
fore data  from  Ref.  11  were  used. 


TEMPERATURE*** 


Figure  1.  Pressure-temperature  regim*  for  the  calculation 
of  air  properties. 
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Figure  2.  Pressure-temperature  regime  for  the  calculation  of 
exhaust  gas  properties. 


The  real  gas  effects  considered  were  intermolecular  forces, 
chemical  dissociation,  vibrationally  frozen  flow,  and  ionization. 

The  effects  of  intermolecular  forces  on  the  constituents  of  air  are 
important  and  are  developed  in  more  detail  in  Section  2 as  well  as 
the  effects  of  chemical  dissociation  in  the  exhaust  gas.  Chemical 
dissociation  of  air  was  considered,  but  the  upper  temperature  limit 
was  below  the  value  where  dissociation  becomes  important.  A series 
of  calculations  were  performed  to  determine  whether  vibrationally 
frozen  flow  should  be  included  in  the  thermodynamic  model.  The 
results  (Ref.  14)  indicated  that  conditions  would  be  highly  unlikely 
to  occur  in  a gas  turbine  engine  where  these  types  of  losses  would  be 
considered  significant.  Ionization  of  air  and  exhaust  gas  was  con- 
sidered, but  since  the  maximum  exhaust  gas  temperature  (4000°R, 
Fig.  2)  is  below  the  temperature  wrhere  ionization  becomes  signifi- 
cant (Ref.  15),  no  corrections  were  included. 
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A number  of  methods  for  computing  thermodynamic  properties  of 
• exhaust  gases  are  available  for  computer  use.  They  fall  into  two 

general  categories;  one  contains  tabulated  combustion  properties, 
while  the  other  is  a computer  program  that  calculates  the  properties. 
Examples  of  the  former  data  format  include  Keenan  and  Kaye  (Ref.  3), 
Banes,  et  al.  (Ref.  6),  and  Powell,  et  al.  (Ref.  16).  Examples  of  the 
latter  include  Pinkel  and  Turner  (Ref.  17),  Osgerby  and  Rhodes 
(Ref.  18),  and  Gordon  and  McBride  (Ref.  19).  The  use  of  the  tabular 
form  of  pre-calculated  data  using  some  form  of  curve-fitting  was  re- 
jected because  of  the  desire  to  retain  control  over  the  air  and  fuel  com- 
position and  thermodynamic  properties  of  the  constituents  without 
having  to  recalculate  the  tables. 

! 

The  thermodynamic  model  must  be  developed  for  on-line  data  re- 
duction use,  which  means  that  calculation  times  and  computer  memory 
requirements  must  be  minimized,  but  sufficient  accuracy  must  be  re- 
tained in  the  constituent  property  data  so  that  the  errors  in  the  calcu- 
i lated  values  are  not  of  the  same  size  as  the  error  in  the  experiment. 

The  maximum  time  per  calculation  (temperature  and  pressure  known, 
to  calculate  enthalpy,  entropy,  or  specific  heat)  should  be  in  the  range 
of  0.05  sec  and  the  memory  requirements  should  be  in  the  range  of 
10,000  words.  By  assuming  the  engine  measurement  error  to  be 
approximately  one  percent,  the  desired  error  between  the  original 
thermodynamic  data  for  the  constituents  and  the  data  calculated  for  use 
in  the  model  should  not  exceed  0. 1 percent  (2  standard  deviations). 

This  report  describes  a computer  model  developed  to  calculate  the 
thermodynamic  properties  of  the  working  fluids  of  a turbine  engine. 

The  model  calculates  the  thermodynamic  properties  of  air  and  the 
exhaust  gas  from  the  reaction  of  air  with  a hydrocarbon  fuel.  Composi- 
tion of  the  air  and  fuel  used  with  this  model  is  discussed.  The  necessary 
thermodynamic  properties  of  the  constituents  are  presented  as  curve  fits 
of  the  data  from  Refs.  10  and  11.  The  logic  involved  in  the  model  is  dis- 
cussed as  well  as  the  method  of  operation.  Calculated  thermodynamic 
data  from  this  model  are  compared  with  various  data  currently  in  use. 

A listing  of  the  model  is  included. 


2.0  DEVELOPMENT  OF  THE  GENERAL  COMPUTER  MODEL 


The  computer  model  calculates  the  thermodynamic  properties 
(enthalpy,  entropy,  specific  heat,  and  Gibbs  free  energy)  for  the  con- 
stituents of  the  working  fluids  of  a turbine  engine,  air  and  exhaust  gas. 
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In  addition,  the  model  calculates  the  enthalpy,  entropy,  specific  heat, 
sonic  velocity,  and  ratio  of  specific  heats  for  air  and  exhaust  gas. 

The  model  will  determine  the  composition  and  molecular  weight  of 
the  exhaust  gas,  but  these  two  quantities  are  fixed  for  air. 

The  model  is  divided  into  two  sections,  the  calculation  of  the 
air  properties  and  the  exhaust  gas  properties.  The  calculation  of 
the  air  properties  includes  a correction  for  intermolecular  forces 
of  the  constituents.  The  calculation  of  the  exhaust  gas  properties 
includes  the  effects  of  chemical  dissociation  at  high  temperatures. 

A schematic  of  the  calculation  process  for  the  model  is  shown  in 
Fig.  3. 
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Figure  3.  Schematic  of  computer  model. 


2.1  CALCULATION  OF  THE  THERMODYNAMIC  PROPERTIES 
OF  DRY  AIR 

2.1.1  Composition  and  Molecular  Weight  of  Dry  Air 

The  composition  of  air  js  based  on  the  U.S.  Standard  Atmosphere 
of  Ref.  9.  This  composition  was  selected  because  of  the  detailed  re- 
search that  went  into  its  original  development  and  the  ongoing  work 
that  can  be  used  to  update  the  composition  if  this  becomes  necessary. 
This  particular  source  is  also  widely  accepted  in  the  turbine  engine 
and  related  industries. 
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The  composition  of  Ref.  9 contains  seventeen  constituents  (see 
Table  1).  To  lessen  storage  requirements  and  computation  time,  the 
number  of  trace  constituents  was  reduced.  The  criteria  for  reducing 
the  number  of  constituents  are  listed  below: 

1.  The  molecular  weight  of  the  final  composition  would 
not  deviate  more  than  ±0.  010  percent  from  the  value 
in  Ref.  9 (28. 9644  lbm/lbm-mole). 

2.  If  the  number  of  constituents  necessary  to  meet  item  1 
is  equal  to  or  greater  than  three,  the  monatomic  trace 
elements  will  be  included  with  a Noble  gas  and  the  re- 
maining trace  constituents  included  with  one  of  the 
other  molecular  constituents. 


Using  a four  constituent  composition  (N2,  O2,  A,  and  CO2)  gives 
a molecular  weight  of  29.9651  if  the  monatomic  elements  are  included 
with  argon  and  the  remainder  with  carbon  dioxide.  This  molecular 
weight  exceeds  the  0. 010 -percent  criteria  and  was  rejected.  A five 
constituent  composition  (N2,  O2,  A,  CO2,  and  Ne)  gives  a molecular 
weight  of  28.9646  lbm/lbm-mole  if  the  remaining  monatomic  con- 
stituents are  included  with  neon,  and  the  other  constituents  are  in- 
cluded with  the  carbon  dioxide.  A summary  of  the  composition  is  shown 
below: 


Constituent 

Mole 

Fraction,  n^ 

Constituent 

Molecular  Weight,  MWj 

ni,  MWj 

Nitrogen 

0.78084 

28.0134 

21.8740 

Oxygen 

0.209476 

31.9988 

6.7030 

Argon 

0.00934 

39.944 

0.3731 

Carbon  Dioxide 
(+  trace 
Constituents) 

0.0003194 

44.00995 

0.0141 

Neon 

(+  monatomic 
trace  constitu- 
ents) 

0.0000246 

20. 183 

0.0005 

1.0000 

28.9646 

The  molecular  weight  of  this  composition  differs  from  that  in  Ref.  9 
by  0.0002  or  approximately  1 part  in  150,000.  Also  the  trace 
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constituents  that  have  to  be  included  with  the  five  constituents  amount 
to  0.00001182  percent. 


The  computer  model  considers  this  composition  and  molecular 
weight  as  fixed.  If  either  the  composition  or  molecular  weight  of  one 
of  the  constituents  were  to  change  significantly,  a simple  card  change 
would  update  the  parameters.  The  inclusion  of  water  vapor  can  also 
be  made  as  a constituent  of  air. 

2.1.2  Equations  for  the  Calculation  of  the  Thermodynamic 
Properties  of  the  Constituents 


The  thermodynamic  properties  of  primary  interest  are  enthalpy, 
entropy,  and  specific  heat  at  constant  pressure.  The  properties  of 
the  constituents  are  taken  from  Refs.  10  and  11.  The  base  tempera- 
ture for  enthalpy  is  taken  as  0°R,  and  the  reference  pressure  for 
entropy  is  one  atmosphere.  The  tabular  data  from  these  sources 
were  put  into  equation  form  for  computer  solution.  The  overall 
range  of  the  data  is  from  300  to  1, 800°R,  and  the  data  were  divided 
into  two  ranges:  from  300  to  900°R  and  from  900  to  1,  800°R.  The 
type  equation  chosen  for  the  data  and  its  degree  was  based  on  keeping 
the  error  between  the  original  and  the  calculated  data  less  than  0. 1 
percent  (2  standard  deviations).  Attempts  were  made  to  fit  the  data  to 
simple  polynominal  equations  such  as 


'iT 


Vi- 


and 


f.  f gjT  + • • • + fyr4 

The  data  in  the  low  temperature  range  could  not  be  fitted  to  the  above 
type  equations  and  meet  the  accuracy  requirements;  therefore  a 
special  technique  discussed  in  Ref.  20  was  used.  This  particular 
method  gives  the  exact  value  at  the  tabulated  points.  Since  Ref.  10 
is  tabulated  in  180°R  increments,  additional  data  were  requested 
from  the  source  (Ref.  10)  in  18-deg  increments  to  provide  a smoother 
curve  fit.  The  data  in  the  900  to  1, 800°R  were  fitted  using  the 
simple  polynomial  expressions  given  below: 
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Entropy 


U.  - a.,  + b,,T  * c.J2  , d;  T3  + e T4 


1 'J  ‘J  >J 


'J  'J 


(1) 


si  « aik  4 - cikT2  ♦ dikT3  + eil(T4 


(2) 


Specific  Heat 


Si  ' "in  * biJ  » «i,T2  * S'!'3  * *i„T4 


(3) 


The  polynomial  curve  fit  equations  were  generally  third  degree  or 
less.  The  polynomial  equations  for  the  constituents  of  air  will  be 
found  in  Table  2 (900  to  1, 800°R  range).  The  fit  of  the  original 
JANAF  data  for  the  major  constituents  gave  a maximum  error  (2  stand- 
ard deviations)  of  0.  1 percent  when  the  data  are  recalculated  for  use  in 
the  computer  model.  Special  precautions  were  taken  to  prevent  a step 
from  occurring  in  the  properties  at  the  transition  region  at  900°R. 

2.1.3  Calculation  of  the  Properties  of  Air 

The  calculation  of  the  enthalpy,  entropy,  and  specific  heat  (Cp) 
is  made  by  summing  the  properties  of  the  individual  constituents 
based  on  their  mole  fractions  in  air.  The  equations  used  are: 

Enthalpy: 


Entropy: 


"a, 


I n j H j 
i=i 


Specific  Heat: 


(4) 


(5) 


(6) 
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The  mole  fractions  of  the  individual  components  of  air  (ni)  and  the 
properties  of  these  individual  components  are  discussed  in  Section  2.1.1 
and  2.1.2,  respectively. 

2.1.4  Correlation  for  Intermolecular  Forces 

Intermolecular  forces  produce  a significant  effect  on  thermo- 
dynamic properties  of  gases  at  higher  pressures  and  lower  tempera- 
tures. This  effect  can  be  seen  in  Fig.  4 where  a nondimensionalized 
form  of  the  specific  heat  at  constant  pressure  (cp/R)  of  air  (from 
Ref.  4)  is  compared  as  a thermally  perfect  gas  and  real  gas.  The  per- 
fect gas  property  is  calculated  from  spectroscopic  data,  and  a virial 
correction  is  applied  to  this  for  the  real  gas  property.  In  Fig.  4,  the 
influence  of  pressure  can  be  seen  as  the  value  of  Cp/R  increases  as 
the  pressure  increases  at  the  lower  temperatures.  To  ensure  that 
the  calculated  values, of  enthalpy,  entropy,  and  specific  heat  for  air 
in  this  model  (Eqs.  (>4),  (5),  and  (6))  include  the  effect  of  intermolecu- 
lar forces,  a correction  is  added  to  the  calculation  of  these  properties 


Figure  4.  Effect  of  pressure  on  a nondimensionalized  specific  heat 
(Cp/R)  for  air  at  low  temperatures. 
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for  each  of  the  constituents  in  air.  The  constituent  corrections  are 
then  added  to  the  properties  calculated  by  Eqs.  (1),  (2),  and  (3). 
These  corrected  constituents  are  then  summed  according  to  Eqs.  (4), 
(5),  and  (6)  to  get  the  properties  of  air  at  the  particular  temperature 
and  pressure. 


The  equations  for  calculating  the  correction  are  taken  from 
Ref,  21  and  are  given  below: 


Enthalpy: 


correct  ion  • 


Entropy: 


correction. 


Specific  Heat: 


(7) 


(8) 


correction 


(t2  d2B  1 f L T d»\2  r dC;  T2  d2cil) 


(9) 


These  corrections  are  then  added  to  the  value  of  the  thermodynamic 
property  calculated  according  to  Eqs.  (1),  (2),  and  (3)  as  given  below: 


Enthalpy: 


Hj  - Hj  + Hcorreclion, 


(10) 


Entropy: 

Sj  = Sj  + Sc0rrectj0n^  (11) 

Specific  Heat: 


correction; 


(12) 
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The  (virial)  coefficients  and  their  derivatives  for  use  in  Eqs.  (7), 
(8),  and  (9)  were  taken  from  data  in  Ref.  22.  These  data  were  fitted 
to  polynomial  equations  for  ease  of  calculations.  The  equations  are 
given  in  Tables  3 and  4. 

2.1.5  Summary  of  Air  Properties  and  Equations 

A five-constituent  air  is  defined  and  shown  below: 


Constituent 

Nitrogen 

Oxygen 

Argon 

Carbon  Dioxide 
Neon 


Mole  Fraction 

0. 78084 
0. 209476 
0. 00934 
0.0003194 
0. 0000246 


The  molecular  weight  is  28.9646.  The  thermodynamic  properties  of 
the  air  (enthalpy,  entropy,  and  specific  heat  at  constant  pressure) 
are  calculated  from  a summation  of  the  properties  of  the  constituents 
based  on  the  composition  (Eni  Hj,  etc.).  The  constituent  properties 
are  curve  fits  of  the  data  from  Refs.  10  and  11.  The  data  are  fitted 
over  a range  from  300  to  1,800°R  using  the  spline  fit  technique  of 
Ref.  20  for  the  300  to  900°R  range  and  polynomial  curve  fits 


from  900  to  1,800°R.  A correction  for  intermolecular  forces 
(Eqs.  (7),  (8),  and  (9))  is  added  to  thermodynamic  properties  calcu- 
lated from  Refs.  10  and  11.  The  equations  for  the  constituent 
properties  900  to  1,800°R  range  and  virial  coefficients  will  be  found 
in  Tables  2,  3,  and  4. 


2.2  CALCULATION  OF  THE  THERMODYNAMIC  PROPERTIES  OF 
TURBINE  EXHAUST  GAS 

The  model  for  calculating  thermodynamic  properties  of  exhaust 
produced  from  the  combustion  of  a hydrocarbon  fuel  (CnHxn)  and  air 
consists  of: 

1.  Determining  the  exhaust  gas  composition. 
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2.  Calculating  the  thermodynamic  properties  of  the 
constituents  of  the  exhaust,  and 

3,  Summing  the  constituent  properties  based  on  the 
mole  fraction  of  constituent  present  in  the  exhaust 
gas. 

The  primary  thermodynamic  properties  calculated  are  enthalpy, 
entropy,  specific  heat  at  constant  pressure,  and  molecular  weight. 

In  addition,  properties  such  as  the  ratio  of  specific  heats  and  the 
sonic  velocity  are  available  as  options. 

’ I 

The  calculations  can  be  made  over  a temperature  range  from 
600  to  4,000°R,  a pressure  range  from  0. 1 to  600  psia,  and  fuel- 
to-air  ratios  from  zero  to  stoichiometric  (Fig.  2).  The  method 
of  calculation  is  dependent  upon  whether  the  effects  of  chemical 
dissocation  are  included.  The  criteria  for  determining  whether 
the  effects  of  chemical  dissociation  are  included  are  based  on  the 
dissociation  of  the  oxygen  molecule  and  are  similar  to  the  criteria 
of  Ref.  23.  If  the  following  equation  is  true 

-i 

T - T|  (13)  ! 

* 

where  j 


i,  = 3000  . 182  lop  Pa 


(14) 


then  the  effects  of  chemical  dissociation  of  the  exhaust  gas  are  not 
considered.  If 


1 1 1 (15) 

the  effects  of  chemical  dissociation  are  considered.  Equation  (14)  is 
similar  to  the  criteria  of  Ref.  24  and  is  based  on  the  amount  of  dissoci- 
ation in  the  oxygen  molecule.  When  T > Tj,  the  mole  fraction  of  the 
dissociated  oxygen  molecule  will  be  0.00500  or  greater.  Figure  5 
shows  the  area  of  interest  for  the  exhaust  gas  calculations  with  the 
division  between  the  area  where  the  effects  of  chemical  dissociation 
are  considered  and  not  considered.  This  divismr.  has  the  benefit  of 
removing  most  turbine  calculations  from  the  area  where  the  effects 
of  chemical  dissociation  are  included. 
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TEMPERATURE  * *R 

Figure  5.  Pressure-temperature  regime  for  the  calculation  of 
exhaust  gas  properties  (considering  dissociation). 

The  computer  model  for  calculating  the  thermodynamic  proper- 
ties uses  the  air  composition  given  in  Section  2. 1 and  a hydrocarbon 
fuel  composition  of  The  fuel  grade  can  be  varied  by  suitable 

choice  of  x"  and  by  changing  the  hydrogen-to-carbon  in  the  model. 
The  thermodynamic  properties  of  all  constituents  for  air  and  exhaust 
gas  are  kept  accessible  in  order  that  properties  may  be  changed  as 
new  data  become  available. 

2.2.1  Thermodynamic  Properties  of  an  Exhaust  Gas  with  the 
Effects  of  Chemical  Dissociation  not  Considered,  T < T| 

The  thermodynamic  properties  for  an  exhaust  gas  that  is  com- 
pletely reacted  with  no  effects  of  chemical  dissociation  considered 
are  calculated  for  T < Tj . The  assumed  reaction  and  constituents 
involved  are  shown  below: 

a On  II  xn  4-  b N.,  ♦ c Oo  * d A t t CO  > * f Ne  -* 

pC.Oo  * j IloO  * I.  No  4 kOo  f dA  * f Ne  (16) 
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1 


! 
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The  reaction  is  assumed  to  take  place  between: 


600 


T 


(17) 


The  thermodynamic  properties  of  the  constituents  in  Eq.  (16)  will 
be  found  in  Table  5.  These  data  are  curve  fits  of  the  tabular  data 
found  in  Refs.  10  and  11  and  cover  the  temperature  range  in  Eq.  (17). 
The  method  of  making  the  exhaust  gas  calculations  is  basically  the 
same  as  Refs.  17  and  25.  The  values  of  the  thermodynamic  proper- 


ties h,  s,  Cp,  and  R of  the  exhaust  gas  are  computed  on  the  basis  of 


the  following  assumptions: 


1. 


The  composition  of  the  exhaust  gas  is  based  on  the  con- 
dition that  the  fuel  is  completely  converted  to  CO2  and 
H2O.  From  this  follows 


The  composition  of  the  exhaust  gas  does  not 
change  in  going  through  a thermodynamic 
process. 


b. 


The  amount  of  unburned  hydrocarbons  in  the 
exhaust  is  negligible. 


2. 

3. 


The  internal  energy  states  of  each  component  gas  are  in 
equilibrium. 


The  exhaust  gas  behaves  as  a perfect  mixture  of  perfect 
gases.  Intermolecular  forces  are  not  considered. 


As  a result  of  assumption  3,  the  thermodynamic  properties  such  as 
enthalpy  of  a mixed  gas  is  the  sum  of  the  enthalpy  of  each  component 
multiplied  by  the  ratio  of  the  mass  of  that  component  to  the  total  mass 
of  the  mixture.  For  up  to  stoichiometric,  the  enthalpy  of  the  exhaust 
gas  is: 


hr  - 


(iTi)[kAi»  * 


(18) 


A*  = 


HHoO  “ 2HOo 


2.016 


IT 


11  CO.,  - «0„ 


12.011 


(19) 

(20) 


21 


-./A 
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The  term  [( Am  + B )/(l  + m)]  in  Eq.  (18)  accounts  for  the  difference 
between  the  enthalpy  of  the  water  and  carbon  dioxide  in  the  reacted 
mixture  and  the  enthalpy  of  the  oxygen  removed  from  the  air  (hair)  by 
their  formation.  The  term  (1/(1  + f)]  expresses  the  units  of  mass  of 
the  enthalpy  in  terms  of  lbm  of  exhaust  gas.  A similar  relation 
between  the  properties  of  the  mixture  and  those  of  the  constituent 
gases  holds  with  regard  to  s,  Cp,  and  R. 

The  molecular  weight  of  the  exhaust  gas  is  determined  by: 


2.2.2  Thermodynamic  Properties  of  Turbine  Exhaust  Gas  When  the 
Effects  of  Chemical  Dissociation  are  Considered,  T > T| 

The  calculation  of  thermodynamic  properties  at  temperatures 
greater  than  Tj  include  the  effects  of  chemical  dissociation 
(see  Fig.  5).  In  Section  2.2.1,  the  exhaust  gas  composition  was  de- 
fined as  containing  CO2,  H2O,  and  constituents  of  air.  When  chemi- 
cal dissociation  becomes  important,  the  exhaust  gas  will  contain 
numerous  trace  constituents.  In  addition  to  those  shown  in  Eq.  (16), 
the  following  are  assumed  as  possible  constituents:  CO,  H2,  O,  OH, 
H,  and  NO.  Other  constituents  might  be  included,  but  the  criteria 
were  to  eliminate  any  constituent  whose  mole  fraction  was  less  than 
10 over  the  envelope  shown  in  Fig.  2. 

The  calculation  of  a thermodynamic  property  such  as  enthalpy, 
entropy,  or  specific  heat  is  made  in  two  steps  when  the  known  value 
of  temperature  "T"  is  greater  than  Tj.  The  procedure  is  to  calcu- 
late the  sensible  enthalpy  at  Tj  by  the  method  of  Section  2.2.1,  then 
to  calculate  the  change  in  chemical  plus  sensible  enthalpy  going  from 
Tj  to  T considering  the  exhaust  gas  to  have  a variable  composition. 
The  enthalpy  change  going  from  Tj  to  T is  then  added  to  the  enthalpy 
at  Tj  to  give  the  enthalpy  at  T. 

The  calculation  of  the  thermodynamic  property  in  the  region 
from  Tj  to  T where  the  composition  is  variable  is  done  by  the  method 
of  Ref.  18.  The  basic  calculation  method  involves  calculating  the 
exhaust  gas  composition  by  the  minimization  of  free  energy.  With 
the  composition  at  Tj  and  T,  the  enthalpy  is  then  calculated  from  the 
Eq.  (22): 


22 
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’(T.P) 


i=,  nH*(Tj,P) 
MW(T.P) 


The  enthalpy  is  also  calculated  at  (Ti„P): 


niHi(T,,P) 

'(T«'P)  = MVi,,p) 


(22) 


(23) 


The  change  in  enthalpy  is  then  the  difference  between  Eqs.  (23)  and 
(22)  and  is  the  sensible  component  of  the  enthalpy  difference: 


Ah(T,P)  = h(T,P)  “ h(T,.P) 

(^ 

The  enthalpy  calculated  for  the  frozen  composition  at  Tj  by  the 
method  of  Section  2.2.1  (Eq.  (18))  is  added  to  the  sensible  enthalpy 
of  Eq.  (24): 


H{TP)  = h(T^  + 


Ah 


(T.P) 


(25) 


The  calculated  value  of  enthalpy  is  always  a positive  value.  Using  the 
method  of  Eq.  (25)  eliminates  the  problem  of  a discontinuity  at  Tj  and 
there  will  be  a smooth  transition  between  the  two  calculation  steps. 
The  calculation  for  entropy  and  specific  heat  at  constant  pressure  is 
made  in  the  same  manner.  The  calculation  of  the  molecular  weight 
is  based  on  the  composition  calculated  at  T and  P.  The  molecular 
weights  summed  as  shown  below: 


\W 


(TP) 


k 

I n;MWj 
i«:  ' ' 


(26) 


The  thermodynamic  properties  of  the  various  constituents  needed 
for  the  model  will  be  found  in  Table  5.  These  properties  are  poly- 
nomial curve  fits  of  the  data  from  Refs.  10  and  11.  The  equations 
are  generally  third  degree  or  less. 


2.2.3  Summary  of  Equations  for  the  Calculation  of 
Exhaust  Gas  Properties 

A computer  model  was  presented  that  calculates  the  thermo- 
dynamic properties  of  the  exhaust  gas  formed  from  the  reaction  of 
air  and  hydrocarbon  fuel  (CnHxn).  The  model  covers  the  tempera- 
ture range  from  600  to  4,00(ER,  the  pressure  range  from  0.1  to 
600  psia,  and  a fuel-to-air  ratio  from  0 to  approximately  stoichio- 
metric. The  model  regime  is  divided  into  two  calculation  regions: 
one  region  where  the  exhaust  gas  is  completely  reacted  into  H2O 
and  CO2  with  no  effects  of  chemical  dissociation  considered,  and 
the  region  where  the  effects  of  dissociation  are  considered.  The 
regions  are  divided  by  the  equation  (Eq.  (14)): 

T,  3000  , 182  log  Pa 


which  is  based  on  the  amount  of  molecular  oxygen  that  dissociates 
in  atomic  oxygen. 

The  thermodynamic  properties  (h,  s,  Cp)  for  the  exhaust  with 
no  effects  of  chemical  dissociation  considered  are  based  on  the 
method  of  Refs.  17  and  25.  The  equation  (Eq.  (18))  for  calculating 
the  enthalpy  is: 

- (t77)[-.  ' '(^1 

and  the  calculation  for  entropy,  specific  heat,  and  the  gas  constant 
follows  the  same  format.  The  method  for  calculating  the  same 
properties  in  the  region  where  the  effects  of  chemical  dissociation 
are  considered  uses  a combination  of  the  method  given  above  and 
the  method  of  Ref.  18.  The  property,  enthalpy  for  instance,  is 
calculated  as  a nondissociated  gas  at  temperature  Tj.  The  compo- 
sition of  the  exhaust  gas  is  then  calculated  at  T and  Tj  (and  the 
known  pressure),  and  the  enthalpy  is  calculated  at  the  two  tempera- 
tures based  on  the  mole  fraction  of  the  various  constituents.  The 
difference  in  enthalpy  at  the  two  temperatures  (T  and  Tj)  is  added 
to  the  nondissociated  enthalpy  calculated  at  Tj  to  give  the  enthalpy 
at  the  known  temperature  and  pressure.  The  calculation  of  the 
entropy  and  specific  heat  follows  the  same  forma';.  The  calculation 
of  the  molecular  weight  is  based  on  a summation  of  the  molecular 
weights  of  the  constituents  weighted  by  the  mole  fractions  calculated 
at  the  known  temperature  and  pressure. 
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The  thermodynamic  properties  of  the  constituents  are  presented  as 
curve  fits  o"  the  tabular  data  found  in  Refs.  10  and  11. 


3.0  OPERATION  OF  THE  COMPUTER  MODEL 


The  computer  model  calculates  the  thermodynamic  properties  of 
air  and  of  the  products  of  combustion  of  air  and  a hydrocarbon  fuel 
(Cn  Hxn).  The  input  parameters  for  the  model  are  temperature  (°R), 
pressure  (psia),  fuel  heating  value  (Btu/lbm),  fuel-to-air  ratio  (lbm/lbm), 
exhaust  gas  enthalpy  (Btu/lbm),  and  XKON  (dimensionless  control  param- 
eter). The  normal  calculation  procedure  is  to  solve  for  enthalpy,  entropy,  J 
spe-  ific  heat,  and  molecular  weight  when  the  temperature,  pressure,  and 
fuel-to-air  ratio  are  known.  The  model  will  calculate  the  equilibrium 
exhaust  gas  temperature  if  the  pressure,  fuel-to-air  ratio,  and  the 
exhaust  gas  enthalpy  are  known.  The  same  calculation  can  be  made  if 
the  fuel  heating  value  is  known  instead  of  the  exhaust  gas  enthalpy.  The 
particular  calculation  is  determined  by  the  computer  through  the  use  of 
dummy  input  parameters.  The  general  operation  of  the  computer  model 
and  the  specific  model  inputs  will  be  described  below.  A general  flow 
diagram  of  the  program  is  shown  in  Fig.  3,  and  a listing  will  be  found 
in  Appendix  B. 

The  model  is  divided  into  two  sections,  one  calculates  the  air  prop- 
erties and  the  other  calculates  the  exhaust  gas  properties.  The  model 
initially  decides  which  section  is  required  from  the  value  of  the  fuel-to- 
air  ratio  (f).  If  the  value  is  zero,  then  the  calculation  of  thermodynamic 
properties  is  made  for  air;  if  a positive  value  is  presented,  then 
properties  for  exhaust  gas  are  calculated  (if  by  accident  a negative 
value  is  input,  the  program  prints  an  error  message  and  the  value  of 
"f").  For  the  calculation  of  the  thermodynamic  properties  of  air,  the 
required  input  data  are  the  known  value  of  air  temperature  and  pres- 
sure and  zero  for  the  fuel  heating  value,  fuel-to-air  ratio,  enthalpy, 
and  XKON.  With  the  value  of  "f"  equal  to  zero,  the  program  begins 
the  calculation  of  air  properties  in  subroutine  THERMO  where  the 
molar  composition  of  air  is  located  as  well  as  the  molecular  weights 
for  the  constituents  of  air  and  the  fixed  molecular  weight  of  air.  The 
thermodynamic  properties  of  air  are  calculated  in  subroutine  GAS,  and 
the  corrections  for  intermolecular  forces  are  calculated  in  subroutines 
HMIXP,  SMIXP,  and  CPMIXP.  If  the  air  temperature  is  between  300 
and  900°R,  the  constituent  properties  of  air  are  calculated  in  subroutine 
HSCP  which  is  the  special  curve  fitting  routine  of  Ref.  20,  and  if  the 
temperature  is  between  900  and  1800°R,  the  constituent  properties  are 
calculated  from  the  polynomial  curve -fits  in  subroutine  COEFF. 


25  i 

I 





AEOC-TR  76-15 


The  calculation  of  the  thermodynamic  properties  of  the  exhaust 
gas  follows  a path  similar  to  air.  Typical  required  inputs  are  the 
known  temperature,  pressure,  and  fuel-to-air  ratio  and  zero  for  the 
fuel  heating  value,  enthalpy,  and  XKON.  In  this  case,  the  positive  "f" 
indicates  that  an  exhaust  gas  calculation  will  be  made  and  then  a check 
is  made  to  see  if  the  temperature  and  pressure  are  in  the  area  where 
the  effects  of  chemical  dissociation  are  included  (Eq.  (14)  or  Fig.  5). 
If  the  properties  are  in  the  nondissociated  region  subroutine  GAS  is 
called  with  temperature,  and  "f. 11  The  thermodynamic  properties  of 
exhaust  gas  are  calculated  in  this  subroutine  using  the  individual  con- 
stituent properties  calculated  from  subroutine  COEFF.  If  the  re- 
quired properties  are  in  the  dissociated  region,  the  temperature, 
pressure,  and  "f"  are  used  in  subroutine  CHEMEQ  to  set  up  the  atom 
concentration  of  the  individual  constituents  in  the  fuel  and  air.  These 
concentrations  and  the  temperature  and  pressure  are  used  in  sub- 
routine PROP  to  solve  for  the  equilibrium  concentration  of  the  exhaust 
gas  species.  Subroutine  PROP  contains  the  convergence  criteria 
(DELL)  which  is  currently  set  for  10“4  for  N2  and  H2O  and  10"6  for 
O2.  This  means  that  the  partial  pressure  of  II2O,  N2,  and  O2  must 
not  change  by  more  than  10“^  or  10"^  between  successive  iterations 
during  the  solving  for  the  specie  concentration.  Generally  a solution 
is  achieved  within  three  iterations.  The  thermodynamic  properties 
such  as  Gibbs  free  energies  are  calculated  in  subroutine  THERM  from 
constituent  properties  found  in  subroutine  COEFF.  When  the  equilib- 
rium concentration  is  established,  then  the  thermodynamic  properties 
like  enthalpy,  entropy,  and  specific  heat  are  determined  in  subroutine 
THERM  from  constituent  properties  found  in  subroutine  COEFF. 


The  model  will  also  determine  the  equilibrium  exhaust  gas  tem- 
perature for  known  values  of  pressure,  enthalpy,  and  fuel-to-air  ratio. 
Zero  is  input  for  the  fuel  heating  value  and  XKON  and  +1 . 0 for  the 
temperature.  The  data  go  through  subroutine  THERMO  into  CHEMEQ 
where  the  atom  concentrations  of  oxygen,  nitrogen,  argon,  neon, 
carbon,  and  hydrogen  are  determined  and  sent  to  subroutine  PROP 
along  witn  the  pressure  and  enthalpy.  Subroutine  PROP  contains  an 
initial  estimate  of  the  temperature  (currently  set  at  3,006°R)  to  be 
used  to  calculate  the  Gibbs  free  energy  and  ultimately  a first  estimate 
of  the  enthalpy.  This  enthalpy  is  compared  with  the  known  value  of 
enthalpy,  and  a correction  is  made  to  the  estimated  temperature  based 
on  the  difference  in  the  enthalpies.  An  iteration  loop  is  set  up,  and 
the  temperature  is  changed  until  the  partial  pressure  of  II2O,  O2,  and 
N2  converge  to  meet  the  requirements  of  DELL  discussed  earlier  and 
the  temperature  change  is  less  than  (0.  l)(cp)  or  approximately  0.03°. 


* 
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The  model  will  also  determine  the  equilibrium  temperature  if  the 
heating  value  of  the  fuel  is  known.  The  input  for  this  case  would  be 
pressure,  fuel  heating  value,  and  fuel-to-air  ratio.  A value  of  -1.0 
is  input  for  the  enthalpy  and  XKON  and  +1.0  for  the  temperature. 

The  negative  values  have  no  physical  meaning,  but  are  used  to  direct 
the  model  to  specific  subroutines.  The  actual  calculation  procedure 
is  the  same  as  for  the  previous  case  where  the  enthalpy  was  known, 
except  that  the  enthalpy  of  the  fuel  must  be  determined  from  the  fuel 
heating  value  in  subroutine  FUELH,  The  convergence  criteria  for 
H2O,  02,  and  N2  are  the  same  as  discussed  earlier. 

A subroutine  (SONV)  has  been  included  with  the  program  to  calcu- 
late the  ratio  of  specific  heats  and  the  sonic  velocity.  This  subroutine 
uses  inputs  (molecular  weight,  fuel-to-air  ratio,  temperature,  and 
the  specific  heat  at  constant  pressure)  from  the  main  program,  but 
the  subroutine  is  not  called  automatically.  Therefore,  if  these  calcu- 
lated values  are  required,  call  statements  must  be  given  in  the  sub- 
routine where  the  values  are  required. 


4.0  PRESENTATION  OF  RESULTS  AND  COMPARISON  OF  DATA 


Various  thermodynamic  properties  were  calculated  using  the 
computer  program  of  this  report  and  were  compared  with  other  data 
currently  in  use  in  the  turbine  engine  industry.  The  compared  data 
consisted  of  enthalpy  and  specific  heat  at  constant  pressure  for  both 
air  and  turbine  exhaust  gas  as  the  working  fluid. 


4.1  COMPARISON  OF  DATA  FOR  AIR 

Enthalpy  is  a primary  thermodynamic  property  used  in  the  analysis 
of  a turbine  engine.  The  enthalpy  of  air  based  on  the  model  described 
was  compared  with  the  published  values  of  enthalpy  in  Refs.  3,  4,  6, 

7,  and  8.  The  enthalpy  was  compared  by  the  percent  deviation  as  a 
function  of  temperature.  The  percent  deviation  was 

(h^j  — h-p)  * 

Percent  Deviation  = — — x 

(hB  “ T j 


calculated  by 
100  (27) 
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The  data  from  the  previously  listed  sources  are  noted  ( >2  while 

the  data  calculated  using  the  method  of  this  report  are  noted  by 
( )j.  Because  of  the  different  base  temperatures  for  enthalpy 

and/or  the  constituent  data  an  artificial  base  was  used,  and  all  data 
were  compared  as  an  enthalpy  difference  between  the  desired  tem- 
perature and  that  of  1,800°R,  The  base  temperature  was  chosen  at 
1,800  R because  air  can  be  treated  as  a thermally  perfect  gas  at  this 
temperature;  also  this  is  the  maximum  temperature  for  the  air  calcu- 
lations. The  data  were  compared  at  two  pressures,  a low  pressure 
of  0.  1 (Fig.  6a)  and  a high  pressure  of  approximately  40  atm 
(Fig.  6b).  Shown  in  Fig.  6a  is  a comparison  of  the  data  from  this 
program  with  three  sets  of  data  where  air  is  treated  as  a thermally 
perfect  gas  (Refs.  3,  4,  and  8)  and  the  enthalpy  is  a function  of  tem- 
perature only.  In  these  references,  there  are  no  corrections  for 
pressure,  and  the  deviation  between  thermally  perfect  gas  behavior 
and  a real  gas  can  be  seen.  Figure  6a  shows  that  air  properties 
as  calculated  by  this  program  approach  perfect  gas  behavior  at 
1, 100°R.  The  data  also  show  that  the  calculated  enthalpy  data  are 
below  that  for  a perfect  gas  above  1, 100°R  and  are  due  to  differences 
in  composition  and  properties.  Figure  6a  shows  a positive  deviation 
for  the  data  comparison  with  Refs.  6 and  7 which  consider  real  gas 
effects.  This  deviation  indicates  that  the  term  (hjggo  ~ h-p^  is 
larger  than  the  same  value  calculated  by  this  program  (hi goo  - h^h 
and,  therefore,  shows  that  the  correction  for  intermolecular  forces 
in  the  two  references  is  larger  than  the  values  used  in  this  program. 
Figure  6b  shows  the  data  as  calculated  by  this  program  for  high 
pressures  (approximately  40  atmosphere)  to  be  very  similar  to  the 
data  in  Refs.  4,  6,  and  7.  This  means  that  the  corrections  for  inter- 
molecular forces  are  approximately  the  same  at  this  pressure.  The 
correction  again  becomes  less  for  this  program  at  1,100°R. 

Figures  7a  and  b compare  the  calculated  value  of  specific  heat  at 
constant  pressure  (Cp)  with  that  from  Refs.  4 and  6.  The  first  refer- 
ence considers  air  as  both  a real  gas  and  a thermally  pe.  feet  gas, 
whereas  the  latter  reference  considers  air  only  as  a real  gas.  The 
deviations  in  both  cases  are  generally  less  than  0.  2 percent  except 
for  40  atm  where  the  deviation  varies  from  -0.2  to  -0.4  percent. 

The  fact  that  the  deviations  are  negative  indicates  that  the  calculated 
values  of  Cp  are  larger,  and  in  this  case,  the  larger  value  is  due  to 
the  difference  in  air  composition  and  constituent  properties  rather 
than  the  corrections.  This  can  be  seen  from  the  fact  that  the  ther- 
mally perfect  gas  Cp  in  Fig.  7a  from  Ref.  4 has  a negative  deviation 
of  the  same  approximate  value  as  the  real  gas  Cp  up  to  a pressure  of 
10  atm. 
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a.  Data  compared  with  that  from  Ref.  4 
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b.  Data  compared  with  that  from  Ref.  6 
Figure  7.  Percent  deviation  of  the  calculated  value  of  Cp  of  air  at  several 
pressures  and  compared  with  data  from  Refs.  4 and  6. 
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A comparison  of  the  deviation  at  360°R  and  approximately  40  atm 
in  Figs.  7a  and  b shows  a large  difference  between  the  data  of  Refs.  4 
and  6.  This  spread  indicates  the  differences  that  can  occur  in  various 
data  sources  at  low  temperatures  and  high  pressures.  In  this  area, 
the  correction  data  used  in  this  computer  model  places  the  calculated 
data  in  between  these  two  sources  and  gives  additional  assurance  that 
the  correction  terms  are  reasonable. 


4.2  COMPARISON  OF  DATA  FOR  EXHAUST  CALCULATIONS 
4.2.1  Comparison  of  Data  with  Those  in  the  Published  Literature 

The  thermodynamic  data  calculated  by  this  portion  of  the  pro- 
gram were  compared  with  similar  data  available  in  the  published 
literature.  The  data  comparison  was  made  with  enthalpy  since  it  is  a 
parameter  of  major  importance  in  the  data  analysis  of  turbine  engines. 
The  comparison  was  presented  in  terms  of  percent  deviation  as  calcu- 
lated by  Eq.  (28). 


(hT-  hg)^  - (hT-  hg)^ 

Percent  Deviation  = - — x 100 

hT  “ hB\  (28) 

The  data  were  calculated  in  terms  of  enthalpy  differences  be- 
tween the  known  temperature  and  the  base  temperature.  The  base 
temperature  was  600°R  when  data  from  Ref.  15  were  being  compared 
and  540°R  when  data  from  Ref.  6 were  used.  To  compare  the  data, 
fuel  was  assumed  CnH2n  to  be  consistent  with  the  data  in  the  references. 
This  meant  that  the  hydrogen-to -carbon  atom  ratio  (M)  in  subroutine 
THERMO  was  set  at  0. 1678  (2  x 1,008/12.011)  and  that  the  number  of 
carbon  atoms  (XN)  was  set  equal  to  1.0  and  the  number  of  hydrogen 
atoms  (XM)  was  set  equal  to  2.0  in  subroutines  CHEMEQ  and  MFNOD. 
Also  in  subroutine  CHEMEQ,  the  stoichiometric  ratio  was  defined  as 
0,067623  in  the  equation  for  PHI.  The  input  data  were 

1000°R  < T < 4000° R 
p = 14.696  psia 
f = 0.067623  lbmjnej/'lbmair 
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Figure  8 shows  the  percent  deviation  versus  temperature 
from  the  data  of  Refs.  6 and  16.  The  data  show  good  agreement  from 
2,000  to  2,  800°R.  When  the  data  from  Refs.  6 and  16  include  signifi- 
cant corrections  for  chemical  dissociation  (beginning  at  2,800°R),  the 
deviation  begins  to  increase.  This  increase  continues  until  the  tem- 
perature T = Tj.  This  increase  in  deviation  is  caused  by  this  model 
not  considering  dissociation  at  a lower  temperature  (i.e.  2,800°R). 
The  deviation  between  2,800  and  3,000°R  (or  the  temperature  that 
dissociation  is  considered  in  this  program)  becomes  a constant  bias 
for  the  rest  of  the  data  at  higher  temperatures.  The  presence  of  this 
bias  is  evident  from  Eq.  (25)  where  the  enthalpy  at  a high  temperature 
(above  Tj)  is  the  sum  of  the  enthalpy  at  the  temperature  Ti  plus  the 
difference  in  enthalpy  between  T and  Tj.  This  is  also  born  out  in 
Fig.  8 where  the  deviation  is  shown  to  be  essentially  constant  after 
reaching  the  temperature  Ti  (which  is  3,000°R  for  these  data).  A 
second  data  comparison  is  shown  in  Fig.  9 for  the  same  fuel  compo- 
sition and  pressure,  but  the  fuel-to-air  ratio  is  now  0.  033812.  The 
agreement  with  the  data  from  Ref.  17  is  good,  again  showing  the 
characteristic  increase  in  deviation  up  to  3,000°R  and  then  becoming 
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Figure  8.  Percent  deviation  of  the  enthalpy  of  turbine  exhaust  g8s  as  calculated 
by  the  method  of  this  report  and  Refs.  6 and  16  at  a fuel-to-air 
ratio  of  0.067623. 


essentially  constant.  The  agreement  with  the  data  from  Ref.  3 is 
good  until  the  computer  program  begins  to  consider  dissociation,  and 
the  agreement  becomes  very  poor.  The  large  deviations  with  the  data 
from  Ref.  3 are  caused  by  this  source  not  considering  dissociation, 
and  therefore,  the  enthalpy  difference  will  be  smaller  in  Eq.  (28)  and 
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leads  to  a negative  deviation  that  gets  progressively  larger  as  the 
temperature  increases. 


Figure  9.  Percent  deviation  of  the  enthalpy  of  turbine  exhaust  gas  calculated 
by  the  method  of  this  report  and  Refs.  3 and  17  at  a fuel-to-air 
ratio  of  0.033812. 

4.2.2  A Comparison  of  the  Enthalpy  of  Two  Hydrocarbon  Fuel 
Compositions 

A common  hydrocarbon  fuel  grade  in  use  in  the  gas  turbine  indus- 
try today  is  JP-4.  Because  it  is  a mixture  of  petroleum  distillates, 
the  chemical  composition  is  hard  to  define.  Generally,  the  average 
composition  is  around  CnHjt  g5n  to  CnH2  Qn*  To  determine  the  order 
of  difference  in  enthalpy  that  would  occur  because  of  a change  in 
hydrogen -to -carbon  atom  ratio  in  the  fuel,  a series  of  calculations 
were  made  with  each  composition.  The  required  computer  param- 
eters for  CnH2t0n  were  noted  earlier,  and  those  for  the  composition 
CnHl.95n  are 


m = 0.1636  (in  subroutine  THERMO) 

XN  = 1.0  (in  subroutines  CHEMEQ  and  MFNOD) 

XM  = 1.95  (in  subroutines  CHEMEQ  and  MFNOD) 
Stoichiometric  ratio  = 0.06795  (in  subroutine  CHEMEQ) 
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| The  calculations  were  made  at  two  pressures  14.696  and  1.4696 

psia  and  at  the  stoichiometric  fuel-to-air  ratio.  The  results  are  pre- 
sented in  Fig.  10.  The  data  show  the  deviation  to  be  approximately 
constant  at  0.12  percent.  There  is  a slight  additional  deviation 
caused  by  the  change  in  pressure  from  14.  696  to  1. 4696  psia  when 
the  effects  of  dissociation  are  included,  but  this  change  is  small 
when  compared  with  the  deviation  caused  by  the  fuel  composition. 


Figure  10.  Percent  deviation  of  the  enthalpy  as  a function  of  temperature 
using  a fuel  composition  of  CH2  and  CH1.95. 

4.2.3  Effect  of  the  Thermodynamic  Model  on  the 
Calculation  of  Turbine  Afterburner  Efficiency 

In  recent  years,  there  has  been  some  variation  in  the  chemical 
model  used  to  determine  turbine  engine  performance  for  the  case  of 
low  efficiencies  in  afterburner  tests.  A comparison  was  made  to  de- 
termine if  other  suggested  nodels  made  a significant  difference  in  the 
calculated  afterburner  efficiencies.  Accordingly,  four  models  were 
investigated  as  follows: 

1.  No  effects  of  chemical  dissociation,  but  with  pre- 
scribed proportions  of  CO, 

2.  No  effects  of  chemical  dissociation  but  with  pre- 
scribed quantities  of  unburned  fuel  (assumed  to  be 
CH2)  in  the  exhaust  gas, 
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3.  Complete  chemical  equilibrium,  and 

4.  Complete  chemical  equilibrium,  but  with  prescribed 
quantities  of  unburned  fuel  in  the  exhaust. 

The  method  of  comparing  the  thermodynamic  models  used  by 
Dr.  I.T.  Osgerby  for  this  study  and  the  data  are  discussed  in  Appendix  A. 
The  results  indicated  that  the  chemical  model  considered  had  an  insignifi- 
cant effect  on  calculated  afterburner  efficiency.  In  addition,  the  data  also 
showed  that  the  calculated  afterburner  efficiency  was  not  significantly 
altered  by  considering  chemical  dissociation  in  the  exhaust  gas  stream  as 
opposed  to  assuming  that  the  fuel  is  completely  converted  to  CO2  and 
H2O  up  to  an  exhaust  gas  temperature  of  3, 200°R. 

4.2.4  Review  of  the  Data  for  Internal  Consistency 

To  determine  the  interval  consistency  of  the  computer  program,  two 
cross-checks  were  made  using  the  thermodynamic  properties,  enthalpy, 
entropy,  and  specific  heat  of  air  and  exhaust  gas.  One  check  of  the  data 
used  the  definition  of  Cp 


where  the  enthalpy  of  air  and  exhaust  gas  were  calculated  over  a 100°R 
internal  and  compared  with  the  specific  heat  calculated  at  the  average 
temperature  as  shown  below. 

/h(T  t 100)  ~ Ht\ 

\ 100  /„  CP(T  + 100  2) 

: LI loo 

CP(T  + 100  2)  (30) 

The  comparison  for  air  and  exhaust  gas  as  a function  of  temperature  is 
shown  in  Fig.  11.  The  data  for  air  were  calculated  at  14.  696  psia. 
Figure  11a  shows  an  approximate  deviation  of  -0.15  percent  over  the 
temperature  range  from  300  to  1,700°R.  This  deviation  compares 
favorably  with  the  maximum  error  of  0. 1 percent  allowed  for  the  con- 
stituent properties.  The  deviation  at  900°R  indicates  the  approximate 
error  caused  by  the  numerical  process  since  the  curve  fit  of  all  the 
constituent  properties  is  exact  at  this  temperature.  The  same  com- 
parison was  made  for  exhaust  gas  at  an  "f"  of  0.067623  and  14.  696  psia 
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in  Fig.  lib.  The  deviation  varies  from  +0.30  to  -0.25  percent  which 
is  essentially  the  same  as  for  the  air.  The  comparison  of  the  exhaust 
gas  was  made  from  2,  500  to  3,900°R.  This  temperature  range  covers 
both  the  dissociated  and  nondissociated  range  of  compositions. 


T«mperotur*  -*R 


a.  Deviation  of  air 


Ttmperolurt  - *R 

b.  Deviation  of  exhaust  gas 

Figure  11.  (Ah/AT)P  compared  with  cp  as  calculated  by 

this  computer  model  to  measure  internal  consistency. 
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A second  comparison  was  made  to  compare  enthalpy  and  entropy 
of  air  and  exhaust  gas  using  the  relationship: 


(31) 


This  relationship  was  used  by  calculating  enthalpy  and  entropy  over  a 
100°R  temperature  interval  and  calculating  the  deviation.  The  results 
are  shown  in  Figs.  12a  and  b for  air  and  exhaust  gas,  respectively. 


Tamptroturt  -*R 


a.  Deviation  of  air 


b.  Deviation  of  exhaust  gas 

Figure  12.  Deviation  between  the  function  1/T  (Ah/AT)P  and  (As/AT) 
as  a function  of  temperature  to  determine  the  internal 
consistency  of  the  computer  program  of  this  report. 
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The  deviation  for  air  is  approximately  +1  percent.  and  the  maximum 
deviation  for  exhaust  gas  is  approximately  +1. 5 percent.  The  devi- 
ation is  believed  to  be  caused  by  calculating  a series  of  constituent 
properties  and  summing  them  rather  than  any  inconsistency  between 
the  data  for  enthalpy  and  entropy.  This  is  based  on  the  calculated 
deviations  shown  in  Fig.  12a  between  300  and  9003R  where  the  calcu- 
lated constituent  properties  are  in  essentially  exact  agreement  with 
Refs.  10  and  11,  yet  the  sum  of  these  properties  for  the  air  mixture 
shows  a deviation  from  +0.75  to  -1.1  percent.  The  data  for  the 
exhaust  gas  show  a similar  deviation,  and  by  removing  the  end  points, 
the  deviations  for  the  air  and  exhaust  gas  are  similar.  The  size  of 
the  deviations  is  probably  inherent  in  the  calculation  process  and  not 
a function  of  any  internal  inconsistencies  in  the  basic  data. 


4.3  COMPUTER  CALCULATION  TIMES 

One  of  the  objectives  of  the  program  was  to  develop  computer 
software  that  would  not  require  excessive  computer  time.  To  deter- 
mine the  length  of  time  required  to  make  the  various  calculations,  a 
time  check  was  incorporated  in  the  program.  The  time  required  to 
make  the  calculations  (T,P,  and  f known  with  h,  s,  and  Cp  to  be  calcu 
lated)  are  listed  below: 

Constituents  Time,  sec 


Air  (with  intermolecular  corrections  0.009 

Exhaust  gas 

nondissociated  0.003 

dissociated  0.028 

These  calculations  were  determined  using  an  IBM  370/155  computer. 

5.0  CONCLUSIONS 

A computer  model  was  developed  to  calculate  the  thermodynamic 
properties  of  air  (300°R  < T < 1800°R,  0,007  atm  < p < 40,  8 atm, 
and  the  exhaust  gas  (600°R  < T < 4000°R,  0.007  atm  < p < 40.8  atm, 
0 < f/ ^stoichiometric  - l.b)  of  a turbine  engine.  The  composition  of 
the  air  was 
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Constituent 


Mole  Fraction 


N2 

0. 78084 

o2 

0.209476 

A 

0.00934 

C02 

0.0003194 

Ne 

0.0000246 

and  the  fuel  composition  was  ^11^.  Several  real  gas  effects  were 
investigated,  and  a correction  for  intermolecular  forces  was  applied 
to  the  constituents  of  air.  A correction  for  chemical  dissociation 
was  applied  to  the  exhaust  gas. 

The  thermodynamic  data  for  the  constituents  of  air  and  exhaust 
gas  were  taken  from  the  JANAF  tables  (Ref.  10)  and  curve  fitted  for 
use  in  the  computer  model.  The  data  as  calculated  for  use  in  the 
model  had  a maximum  error  of  0. 1 percent  (2  standard  deviations) 
when  compared  with  the  JANAF  tabulations.  The  calculated  enthalpy 
data  for  air  were  compared  at  0. 1 and  40  atm  with  NBS  data  (Ref.  4). 
The  maximum  deviations  expressed  as 

(hy  — ^ I800°^N BS  “ (hT  “ h1800°^Model 

— — 7T j— 1 x 100 

'hT  “ h1800'Model 

were  approximately  -0.08  and  -0. 10  percent,  respectively.  The 
calculated  enthalpy  for  a turbine  exhaust  gas  (CnHxn/Air  system  at 
14.  696  psia  and  a fuel-to-air  ratio  of  0.0676)  was  compared  with 
similar  data  from  the  AGARD  tables  (Ref.  6).  The  maximum  devi- 
ation (expressed  in  a similar  manner  to  air)  was  0.4  percent. 

A set  of  calculations  was  performed  to  see  if  the  low  efficiencies 
frequently  experienced  during  afterburner  testing  could  be  traced  to 
the  chemical  model  employed.  The  data  showed  that,  for  all  models 
considered,  the  calculated  efficiencies  were  essentially  the  same. 

Calculation  times  were  determined  for  different  type  data. 
Typical  times  for  a known  T and  P are  listed  below: 
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Constituents 


Time,  sec 


Air  (with  intermoiecular  corrections)  0.009 

Exhaust  gas 

nondissociated  0.003 

dissociated  0.028 

These  calculations  were  determined  using  an  IB]\]  370/155  computer. 
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Table  1.  Air  Composition  from  Ref.  9 


Constituent  Gas  Content*,  percent  by  volume 


1. 

Nitrogen  (Nj) 

78.084 

2. 

Oxygen  (Oj) 

20.9476 

3. 

Argon  (Ar) 

0.934 

4. 

Carbon  dioxide 

(co2) 

0.0314 

5. 

Neon  (Ne) 

0.001818 

6. 

Helium  (He) 

0.000524 

7. 

Krypton  (Kr) 

0.000114 

8. 

Xenon  (Xe) 

0.0000087 

9. 

Hyrdogen  (Hj) 

0.00005 

10. 

Methane  (CH^) 

0.0002 

11. 

Nitrous  oxide 

(n2o) 

0.00005 

12. 

Ozone  (0_) 

Summer:  0 to  0.000007 

»»  * t A . _ A A A A A A O 

Winter:  0 to  0.000002 


! 

t 


13. 

Sulfur  dioxide  (S02) 

0 

to 

0.0001 

14. 

Nitrogen  dioxide  (N02> 

0 

to 

0.000002 

15. 

Ammonia  (NH^) 

0 

to 

trace 

16. 

Carbon  monoxide  (CO) 

0 

to 

trace 

17. 

Iodine  (Ij) 

0 

to 

0.000001 

*Numbers  1 through  11  are  included  in  model; 
numbers  12  through  17  are  assumed  to  be  zero. 
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Table  2.  Equations  for  the  Thermodynamic  Properties  of  the  Constituents  of  Air 

Equations  for  Sensible  Enthalpy  (ETU/lba-aola) 

1.  Teaperature  Range  900*-1800*R 

N2>  B - 272.57872  + 6.28498S8T  ♦ 4.2460662  X 10~4T2 
Oj,  B - -49.9253121  ♦ 6.800S911T  ♦ 4.6097941  X 10_4T2 
A , B - 0.19368  + 4.96821S2T 

COj,  B - -291.6316413  ♦ 6.5377190T  ♦ 28.320187  X 10_<T2  - 0.3858033  X 10-ST3 
fc , H - 0.19368  ♦ 4.96821S2T 

Equations  for  Entropy  (less  - RLnp^  terra) , (BTU/lbm-mole,  °R) 

Temperature  Range  900-1800*R 

N2,  S = 39.251394  + 0.0155587T-0. 5615997  X 10_5T2  + 0.933764  X 10*9T3 

02,  S =■  42.269726  + 1.5770590  X 10-2T-0 . 5390427  X 10*5T2  + 0.8573795  X 10_9T3 

A , S = 31.5883510  + 1.2918283  X 10_2T-5. 3921089  X 10_6T2  + 0.9663917  X 10~9T3 

C02,  S =*  41.78458  + 2.0871316  X 10_2T-0. 6347125  X 10*5T2  + 0.9526361  X 10-9T3 

Ne,  S * 29.5537807  + 1.2917042  X lo'2T-5. 3907410'2T-5. 3 907410  X 10_6T2  + 0.9659557  X 10_9T3 

Equations  for  the  Specific  Heat  at  Constant  Pressure  (BTU/lbn-mole*R) 

T— gsrature  Range  900*-1800*R 

Cp  - 7.23667810  - 1.3015396  X 10'3  T ♦ 1.5759008  X 10"*T2  - 0.3746857  X 10_*T3 
Oj,  Cp  ■ 5.7846285  ♦ 2.2415270  X 10-3  T - 0.4580023  X 10*6T2 
A , Cp  - 4.96815 

C°2«  cp  “ 5.3131429  ♦ 0.8363827  X 10’2  T - 0.3088999  X 10*ST2  ♦ 0.4493104  X 10"*T3 
■ a.  Cp  . 4.96815 
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Table  3.  Equations  for  the  Second  Virial  Coefficient  and  the  Derivatives 
for  the  Constituents  of  Air 


1.  Nitrogen  (Nj) 

5.1364637  X 102  -0.6636122T  -0.3534783  X 10~3T2 
B " 1.0-1.7153299  X 10'2T  -0.9296798  X 10*5I2 

T dB  . -3.4560792  X 102 

dX  1.0-1.1286423  X 10"2T  -0.29627746  X 10_5T2 

X2^»  - 6.7503388  X 102  

dT  1.0-0.8951382  X 10"2T  -0.4815325  X lO'^T2  ♦ 1.3834608  X 10*»T3 

2.  Oxygen  (02) 

4.2857084  X IQ2  -4.7667818  X 10-1T  -1.5050447  X 10~4T2 
* " 1.0  -1.4123010  X 10'2T-0. 5095962  X lO-2^ 

_dB  _-3. 1691728  X 102 

Tdf  1:6-6.9511198  x 1'6-7t-0.  JJT833J  xT6*:>tZ  ♦fiTTCimS  X 16-9T3  “ 

T2d2B  m -3.3796998  X 101-4.9627303T  ♦ 1.0935667  X IQ-2?2 
dT2  1.0-1.9686162  X 10*2T  *1.2893683  X 10“4I2 

-0.5888087  X 10~5X3  ♦0.3089713  X 10~*T4 
-0.2105136  X 10*«I3  ♦ 0.946670  X 10_1°T4 

3.  Argon  (A) 

B - 3.8904656  X 102-0.870535T  ♦0.4772262  X I0~3X2  -1  .4860291  X IQ^T3 
lVo-1 . 4845725  X 10_2T+1 . 1790085  X 10_5I2  ♦ 1.4973286  X 10*10T3 

dB  . -2.7524532  X 102  ♦ 4.4273617  X 10~2T 

1.0-0.9485511  X 10*2T 

X2  d^B  - 3.9858659  X 102-3. 4251675  X 10~2T 
dT2  1.0  ♦ 0.6811836  X 10_2T 

4.  Carbon  Dioxide  (COj) 

5.3152906  X 102-5. 0336843  X 10-1T  ♦0.702646*  X 10~*T* 

1.0- 0.6113584  X 10"2T  -2.2018674  X 10_,T* 

5.3458176T  -0.9807759  X IQ'^T  -7.8788948  X 10~*T2  18.9162693  X 10~12T3 

(1.0-0.6113584  X 10_2T-2. 2018674  X 10"8*2 

11. 8074171T2-125. 9326564  X 10~5T*  -1.4462917  X 10~7T4-47. 4547667  X IQ*1 3T*-50 . 1 1 59356  X 
(1.0-1.104452  X 10"2T-7. 13405040  X 10' 

10*7T* 

59.425695  - 0.27144987T" 

1.0- 1.646348  X 10-2t 
-6.1992870  X 101  *1.8886040  X 10~2t 

1.0-1.7014989  X 10’2T 

1.5852079  X 102  -3.1378493  X 10~2T 

1.0- 1.8918656  X 10_2T 


( 1 * Derivative  taken  fro*  fitted  value  of  ‘A* . 


B - 

(1) 


dB 
dT 

T2  dV3*. 
dT 


5.  Neon  (Ne) 

B * 
dB 

T3x  ‘ 

t2  24  - 

dT2 
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Table  4.  Equations  for  the  Third  Virial  Coefficients  and  the  Derivatives 
for  the  Constituents  of  Air 


1 .  Nitrogen  (N2) 

355.28210  - 4. 7253632T 
C - 1.0  - 4.7446815  X 10"  3 T 

dC(1)M  -30.396628  X 10-1T  

3t  “ 1.0  - 9.4893630  X 10*3  T «■  22.5120025  X 10“6T2 

t2  dfc(1,„  -28.8444637  X 10~3  T2  + 136.8577931  X IQ^T3 

dT2  (1.0  - 9.4893630  X 10"3  T + 22.5120025  X 10~6  T2)2 


2.  Oxygen  (02) 

-4.2345758  X 103 

c " 1.0  - 0.9728086  X 10-3T 

t dc<1>  . 4.1194317  X 10-1T 

3f  1.0  - 0.9728086  X 1<T2T 

2 d2Cll)  , -4.0074185  X 10~3T2 
^ " 1.0-0.9728086  X 10“2T 

3.  Argon  (A) 


C - 5.987.9538  - 19.016416T  ♦ 2.4336141  X 10"2I*  -1.0787834  X 10-5I3 

T . -19 .016416T  ♦ 486.72282  X 10  * T - 3.2363502  X 10  3TJ 

2 d2C^^  —4  2 —5  1 

T2  - 486.72282  X 10  * T - 6.4727004  X 10  3T 

dT* 


4.  Carbon  Dioxide  (COj) 

C - 11517.385  - 13. 3839 8T  ♦ 0.3190838  X 10"2T* 

dC^^  -2—2 

- -13.38398T  ♦ 0.6381676  X 10  2T* 

2 d2C^^  -2  2 

T2  — i . 0.6381676  X 10 

dT2 

5.  Neon  (Ne) 


C - 0.0 

T a?  * 0,0 

t’  si  - ... 

dT2 


(1)  Derivative  taken  from  fitted  value  of  C. 
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Table  5.  Equations  for  the  Thermodynamic  Properties  for  the  Constituents 
of  the  Exhaust  Gas 

Equations  for  Sensible  Enthalpy  (BTU/lbm-moXe) 

Temperature  Range  -900*-1800*R 

«2,  U - 272.57872  ♦ 6.2849888T  + 4.2460662  X 10_4T2 
02,  H - -49.9253121  ♦ 6.8005911T  + 4.6097941  X 10_4T2 
A . H - 0.19368  + 4.9682152T 

C02,  H - -291.6316413  «•  6.5377190T  + 28.320187  X 10-4T2  - 0.3858033  X 10-6T3 
Ne,  H - 0.19368  + 4.9682152T 
H20,  H - 356.6510897  + 6.913576T  ♦ 8.1349063  X 10_4T2 
Temperature  Range  - 1800*R-2700*R 

N2,  H - -220.5134204  + 6.8101302T  2.8505033  X 10"4T2 

02,  H - -630.6330969  + 7.3851082T  + 2.599229  X 10_4T2 
A , H - 0.0615  + 4 .9681 41 8T 

C02,  H - -3294.0760689  ♦ 11.07844T  5.4166925  X 10_4T2 

Ne,  H - -0.0615  ♦ 4 .968141 8T 
HjO,  H - 152.8596437  + 7.1089372T  + 7.6785528  X 10_4T2 
Temperature  Range  -2700*-4000*R 

N2,  H - -1304.1576077  ♦ 7.6123722T  + 1.3657170  X 10_4T2 
02,  H - 1204.9978428  ♦ 7.8839888T  + 1.5887864  X 10_4T2 
A , H * -0.02187  ♦ 4.9681271T 
C02,  H - -5537.3281389  ♦ 12.747913T  ♦ 2.3106225  X 10-4T2 
Ne,  H - -0.02187  ♦ 4.9681271T 
H20,  H - 1858.4935816  ♦ 8.5657282T  ♦ 5.0420885  X 10"4T2 

Equations  for  Chemical  plus  SenaibXe  EnthaXpy  (BTU/Xba-moXe) 

Temperature  Range  1800*-2700*R 

N2,  H - -20.5134204  ♦ 6.8101302T  ♦ 2.8505033  X 10*4T2 
02,  H - -630.6330969  + 7.3851082T  ♦ 2.5992293  X 10“4T2 
A , H —0.0615  _ 4.9682152T 
C02,  H - -175045.8247470  ♦ 13.506027T 
Ne,  H - -0.0615  + 4.9682152T 

HjO,  B - -102632.3476190  ♦ 7.10M046T  + 7.676647  X 10~4T2 
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Table  5.  Continued 

Temperature  Range  - 2700*-4000*R 

N2,  U —1304.1576077  ♦ 7.6123722T  ♦ 1.3657170  X 10-4T2 
02,  K —1240.9978428  ♦ 7.88398881  ♦ 1.5887864  X 10-4T2 
A , U - -0.02187  ♦ 49681271T 

COj,  U - -177175.4038470  ♦ 14.294766T 
Ne,  U - -0. 02187  ♦ 4.9681271T 

HjO,  U - -104643.9857996  ♦ 8.56583731  ♦ 5.0418784  X 10-4I2 
'eaperature  Range:  2000*-4000*R 

H,  H - 9.2937814  X 104  ♦ 4.9679289T 

O,  H - 1.0650436  X 105  ♦ 4.9819292T 

OH,  H - 1.7243891  X 104  ♦ 6.5117725T  ♦ 2.4740577  X 10-4T2 

NO,  H - 3.7810593  X 104  ♦ 7.7703187T  ♦ 1.3914575  X 10~4T2 

CO,  H - 4.9950257  X 104  ♦ 7.49910571  ♦ 1.6384251  X 10"4I2 

Hj,  H - 4.379566  X 102  ♦ 6.2305090T  ♦ 2.7323736  X 10-4I2 

3.  Equations  for  Entropy  (less  - RLnpj,  term),  (BTU/lbm-mole,  *R) 

Temperature  Range  900-1 800*R 

N2,  S - 39.251394  ♦ 0.0155587T-0. 5615997  X 10_5T2  ♦ 0.933764  X 10-9T3 

Oj,  S - 42.269726  ♦ 1.5770590  X 10-2T-0. 5390427  X 10-ST2  ♦ 0.8573795  X 10_9T3 

ft  , S - 31.5883510  ♦ 1.2918283  X 10~2T-5. 3921089  X 10_6T2  ♦ 0.9663917  X 10*9T3 

COj,  8 - 41.78458  ♦ 2.0171316  X 10_2T-0 .6347125  X 10_ST2  ♦ 0.9526361  X 10~9T3 

Ne,  S - 29.5537807  ♦ 1.2917042  X 10_2T-5. 39074 10~2T-5. 39074 10  X 10~6T2  ♦ 0.9659557  X 10_9T3 
H2Ot  S - 37.437166  ♦ 1.1092773  X 10”2T-0. 6382909  X lO”5*2  * 1.0746317  X 10”9T3 

Temperature  Range  - 1800*-2700*R 

Nj,  S - 57.275024-1.8300467  X 10-2  T ♦ 1.7559820  X 10_5T2-0. 581 5832  X 10”*T3 

0Jf  S - 47.225630  ♦ 0.7577402  X 10-2T  - 0.828194  X 10-6T2 

A , S - 36.4464480  _ 4.5571014  X 10”3T  - 0.5119033  X 10~6T2 

COJt  S - 47.518398  _ 1.1512786  X 10_2T-1 .2066643  X 10_6T2 

Ne.  S - 34.4116312  ♦ 4.5570565  X 10-3  T-0.51 18049  X 10_6T2 

HjO,  S - 43.261707  ♦ 0.8244883  X 10*2  T - 0.7782285  X 10_6T2 
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Table  5.  Continued 


Temperature  Range  - 2700*-4000*R 


M2,  S - 46.8564740  * 5.0506990  X 10‘3  T - 0.3700634  X 10-6T2 


02,  S - 50.2049115  5.2780447  X 10-3  T - 0.3852719  X 10-6T2 


A , S - 38.5124276  ♦ 3.0262745  X 10-3  T - 0.22832958  X 10"*T2 


C02,  S - 51.5035705  ♦ 8.44472  X 10'3  T - 0.6170116  X 10"6T2 


He.  S - 34.8197332  ♦ 4.5503522  X 10-3  T - 0.6902733  X lo’S2  ♦ 0.0462520  X 10“9T3 


H20,  S - 45.6963147  ♦ 0.6353274  X 10-2  T - 0.4115987  X 10-6T2 


4.  Equations  for  the  Specific  Heat  at  Constant  Pressure  (BTU/lbm-mole*R) 
Temperature  Range  - 900*-1800*R 


»-3 


-6„.2 


*l2'  Cp  “ 7.23667810  - 1.3015396  X 10  T ♦ 1.5759008  X 10  dT‘  - 0.3746857  X 10  T 


02,  Cp  - 5.7846285  ♦ 2.2415270  X 10-3  T - 0.4580023  X 1o”6T2 


A , Cp  « 4.96815 


C02,  C - 5.3131429  ♦ 0.8363827  X 10-2  T - 0.3088999  X 10-5T2  + 0.4493104  X 10-9T3 


He.  Cp  - 4.96815 


HjO,  Cp  - 7.5482768  ♦ 4.2815697  X 10-4  T ♦ 0.7157717  X 10_6T2  - 1.3495489  X 10~10T3 


Teaperature  Range  1800*-2700*R 


H2,  C - 5.7766562  ♦ 1.5061552  X 10’3  T - 0.8888113  X 10  "T‘  + 1.0954859  X 10  ,UT 


-6„2 


i»-10_3 


A , c « 4.96815 

P 

CO,,  c ■ 6.680396  * 0.5990060  X 10’2  T - 1.7096379  X 10_6T2  + 0.1809910  X 10"9T3 
* P 


Ne,  c - 4.96815 
P 


»-3 


-6„2 


H,0,  C - 7.1564387  ♦ 0.9855705  X 10  T ♦ 0.4640690  X 10  ”T*  - 1.0002320  X 10  T 

1 P 


irature  Range  - 2700*  4000*R 


H2,  Cp  - 6.5727317  ♦ 0.9099083  X 10* 3T-0. 9596494  X 10-?T2 
02,  Cp  - 7.6813079  ♦ 4.4281245  X lo’S-l  .8738208  X 1o”8T3 


A , C ■ 4.96815 

C02'  Cp  * 10.8521173  + 1.6107530  X 10_3T-1 .7125528  X lo'V 


Me,  Cp  - 4.96815 


HjO,  Cp  - 5.9764766  ♦ 2.58*1179  X 10-3  T - 2.3692669  X lO-7!2 
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Table  5.  Concluded 

Equations  for  the  Gibbs  Free  Energy  Function  (BTU/lbm-nole) 

T— p erasure  Kangs  - 2000*1-4000*1 

H,  r*-  5.3150524  X 10*  - 15.6253583  T - 0.9767035  X 10_3  X 10"3T2  ♦ 0.5610564  X 10"7T* 

O,  F*»  6.0670223  X 104  - 21.8518477  T - 0.9782325  X 10_3T2  ♦ 0.5614114  X 10  ?T3 

OH,  F*-  1.0843812  X 104  - 24.8213338T  - 1.4149044  X TO-3!2  + 0.7304718  X 10~7T3 

HO,  F*-  2.2896290  X 104  - 28. 1298227T-1 .5979361  X 10-3  T2  ♦ 0.8680263  X 10'7T3 

. 4 -3  2 -7  3 

CO,  F - 2.5226419  X 10  - 26. 3634922T-1 .5649177  X 10  J T ♦ 0.8457005  X 10  T 

H2,  F*-  1.4415362  X 103  - 17.8897838T  - 1.3549595  X 10_3T2  ♦ 0.6781319  X 10_7T3 

M,,  F*-  1.2087000  X 103  - 25. 6082461T-1 .5325905  X 10~3T2  * 0.8179650  X 10~7T3 

02,  F*-  1.2814893  X 103  - 27.3792772  T-1. 6271943  X 10_3T  ♦ 0.8783605  X 10_7T3 

C02,  F*-  -9.2608408  X 104  - 27.9578811  1-2.5501816  X 10-3T2  + 1.3506145  X 10* ?T3 

H20,  F*-  5.6574736  X 104  - 25.00726001  - 1.8841889  X 10_3T2  + 0.8839907  X 1o"?T3 
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APPENDIX  A 

INFLUENCE  OF  CHEMICAL  MODEL  ON  AFTERBURNER 
EFFICIENCY  CALCULATIONS 


In  recent  years,  there  has  been  some  variations  in  the  chemical 
model  used  to  determine  turbine  engine  performance  for  the  case  of 
low  efficiencies  in  afterburner  tests.  A comparison  of  several  chemi- 
cal models  was  made  to  determine  the  differences  in  the  calculated 
afterburner  efficiency.  Four  models  were  investigated  as  follows: 


1.  No  effects  of  chemical  dissociation,  but  with  pre- 
scribed proportions  of  CO, 

2.  No  effects  of  chemical  dissociation,  but  with  pre- 
scribed quantities  of  unburned  fuel  in  the  exhaust, 

3.  Complete  chemical  equilibrium,  and 

4.  Complete  chemical  equilibrium,  but  with  prescribed 
quantities  of  unburned  fuel  in  the  exhaust. 


Calculations  were  carried  out  with  the  models  modified  to  include 
observed  (Refs.  26  to  29)  high  levels  of  carbon  monoxide  (CO)  and  un- 
burned fuel  (as  CH2  in  the  model).  In  the  limit,  either  of  these  species 
can  represent  the  effects  of  less  than  ideal  combustion  efficiency.  The 
computer  model  was  basically  a real  gas  Rayleigh  line  heat  addition 
model  with  exhaust  conditions  and  composition  defined  by  the  program 
of  Ref.  18. 


Model  (1)  was  chosen  to  provide  the  baseline  calculations.  The 
calculated  inlet  conditions  for  the  afterburner  are  listed  in  Table  A-l. 
It  was  assumed  that  the  inlet  fuel /air  ratio  would  be  increased  to 
0.06  in  the  afterburner.  The  inlet  temperature  and  pressure  were 
assumed  to  be  1,300°R  and  0.  5 atm,  respectively,  and  the  inlet  Mach 
number  was  0.2. 

Calculations  were  carried  out  assuming  that  a prescribed  pro- 
portion of  the  theoretical  CO2  (without  CO  in  the  exhaust)  is  converted 
to  CO.  The  calculated  nozzle  exhaust  total  pressures  and  tempera- 
tures together  with  the  combustion  efficiency  are  shown  in  Table  A-2. 
The  efficiency  was  evaluated  as  follows: 


fin)Ahf 
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where  hQ  denotes  sensible  total  enthalpy,  f denotes  fuel /air  ratio, 

Ahf,  the  lower  heating  value  of  fuel,  and  subscripts  in  and  ex  denote 
afterburner  inlet  and  nozzle  exhaust  conditions,  respectively.  The 
sensible  enthalpy  was  evaluated  by  subtracting  the  static  enthalpy 
(sensible  plus  chemical)  at  298,  15°K  from  the  total  enthalpy,  mixture 
composition  being  the  same  for  both  enthalpy  calculations. 

Having  determined  the  exhaust  conditions  for  model  (1),  the 
sensible  enthalpy  changes  and  efficiency  were  evaluated  with  models 
(2)  through  (4),  assuming  the  same  nozzle  exhaust  total  temperature 
obtained  with  model  (1).  The  chemical  equilibrium  was  modified  to 
include  unburned  fuel  as  a species. 

The  calculated  afterburner  efficiencies  for  the  four  models  are 
shown  in  Table  A-3,  and  the  concentrations  of  CO  and  unburned  fuel 
used  in  the  calculations  are  shown  in  Table  A-4.  The  efficiencies  in 
Table  A-3  show  no  significant  difference  for  a given  set  of  conditions 
regardless  of  the  calculation  model  used.  This  means  that  the  calcu- 
lated enthalpies  for  the  various  models  are  essentially  the  same,  and 
therefore,  the  choice  of  models  is  not  important  over  the  range  of 
these  calculations. 


The  ratio  of  efficiencies  for  the  different  models  is  so  close, 
even  at  the  3,2466R  temperature,  that  it  should  be  possible  to  increase 
the  base  temperature  below  which  assuming  no  dissociation  should  be 
an  adequate  model.  This  could  result  in  a significant  reduction  in 
data  manipulation  time  with  on-line  data  acquisition  systems. 

Table  A-1.  Afterburner  Inlet  Conditions 


Parameter 


Total  Conditions  Static  Conditions 


Pressure 

Temperature 

Specific  Heat 

R«tio  of  Specific  Heats 

Entropy 

Speed  of  Sound 
Mach  Number 
Density 

Molecular  Weight 
Fuel-to-Air  Ratio 


0.5  atm 
1300  0 R 

0.267  Btu/lb°R 
1.346 

1.9625  Btu/lb°R 
1733.3  ft/sec 
0 

0.01525  lb/ft3 
28.956  lb/lb-mole 
0.025  lbfuel/lbair 


0.4868  atm 
1291. 1°R 
0.267  Btu/lb°R 
1.346 

1.9625  Btu/lb° R 
1727.1  ft/sec 
0.2 

0.01495  lb/ft3 
28.956  lb/lb-mole 
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Lower  heat  of  combustion  Ah,  « 18,650  Btu/lb 
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Table  A-3.  Afterburner  Efficiencies 


Po 

ex 

T 

o 

ex 

nl 

n2  & n4 

n3 

rij^ /n3 

Vn3 

0.4875 

2567 

0.639 

0.647 

0.643 

0.994 

1.006 

0.4867 

2644 

0.679 

0.688 

0.683 

0.994 

1.007 

0.486 

2720 

0.719 

0.729 

0.723 

0.995 

1.008 

0.4852 

2797 

0.759 

0.768 

0.763 

0.995 

1.008 

0.4845 

2873 

0.799 

0.808 

0.803 

0.995 

1.006 

0.4838 

2948 

0.839 

0.846 

0.843 

0.995 

1.004 

0.483 

3023 

0.879 

0.886 

0.882 

0.996 

1.004 

0.4823 

3098 

0.919 

0.925 

0.922 

0.997 

1.003 

0.4816 

3172 

0.959 

0.964 

0.962 

0.997 

1.002 

0.4809 

3246 

1.000 

1.000 

1.000 

1.000 

1.000 

Table  A 4 

. CO  and  Fuel  Vapor 

Exhaust  Concentrations 

1 

2 

3 

4 

T 

o 

% CO 

ppmv  CO 

ppmv  CO 

& Fuel 

ex 

2567 

5.14 

1.0 

1.7 

1.56 

2644 

4.58 

2.3 

3.4 

1.38 

2720 

4.02 

4.0 

6.4 

1.2 

2797 

3.46 

10.0 

11.8 

1.02 

2873 

2.89 

19.0 

21.0 

0.85 

2948 

2.32 

32.0 

36.1 

0.68 

3023 

1.74 

50.0 

60.3 

0.51 

3098 

1.17 

82.0 

97.8 

0.34 

3172 

0.584 

130.0 

154.7 

0.17 

3246 

0 

239.2 

239.2 

0 
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COMPUTER  PROGRAM  LISTINGS 
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FORTRAN  IV 


0001 

0002 

0003 

0004 

0005 

0006 
000T 
0008 


0009 

0010 
0011 
0012 
0013 
001 A 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 


G LEVEL  21  MAIN  DATE  ■ 75224  13/41/03 

C MAIN  PROGRAM 

REAL  M(MWA1R 
R£A0(5.10)LP«NP 
10  F ORMA T(12«2X«12) 

ITEM  » 0 
16  HR1TE(LP<4) 

4 FORMAT (1M1) 

WRITE  <LP«5) 

5 FORMAT (1H  .35HTHERMOOYNAM1C  PROPERTIES  - CHECKOUT// 120M  TEMPERATUR 

IE  PRESSURE  HEAT.  VALUE  FUEL/AIR  ENTHALPY  ENTROPY  SPE 
2C1F1C  HEAT  WTM  CASE  /> 

15  CONTINUE 

READ (NP, 20) TEMPI, PRES1. HEATV1, FA1R1. HI, XKONl 
20  FORMAT(6ElO.O) 

TEMP  ■ TEMPI 
IF (TEMP)50»50,18 
18  CONTINUE 
PRES  ■ PRES  1 
HEATV  ■ HEATVl 
FAIR  ■ FA1R1 
H * HI 

XKON  > XKONl 
ITEM  * 1TEM»1 

CALL  THERHOITEMP, PRES, HEATV. FAIR, XKON, H,S. CP. WTM.LP. ITEM) 

WR1 TE (LPf  30 ) TEMP,PRES,H£ATV ,FA1R,H,S  * CP, WTM, 1 T£M 
30  FORMAT ( 1H  ,4F12.4,2F12.5,3X,2F12,5,5X,I4) 

45  CONTINUE 
GO  TO  15 
50  CONTINUE 
STOP 
END 
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FORTRAN  IV  G LEVEL  21 


THERMO  DATE  * 75224  13/41/03 


f 

5 

i 


i 


$ 

rr 


0001 

0002 

0003 

0004 

0005 

0006 
0007 

oooe 

0009 


0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 


SUBROUTINE  THERMO (TR,P, MV, F,XKON,HH,SS»CPP»WTM,LP, ITEM) 

REAL  M.MWA1R 

C0MM0N/CATHY/XN2.X02, XAR, XC02.XNE  ,MWA1R,Y(I4),X(I0) * PAR* PNC 
PA  a P/14.696 

C MOLE  FRACTIONS  OF  CONSTITUENTS  OF  AIR 
XN2  a 0.78084 
X02  » 0.209476 
XAR  ■ 0.00934 
XC02  ■ 0.0003194 
XNE  • 0.0000246 

C MOLECULAR  VEIGMT  OF  CONSTITUENTS 

C 1-N?,?.02,3-AR,4-C02.5-NE,6-H20.7-H,8-0,9-OM,10-N0.11-C0,12-M?,13-C,14-N 

Y(l)  ■ 28.0134 
Y (2)  * 31.9988 
Y ( 3)  a 39.944 
Y (41  a 44.00995 
Y<5>  ■ 20.183 
Y(6)  * 18.016 
Y (71  * 1.00797 
Y(8)  « 16.0 
Y ( 9)  a 17.0074 
YI10)  a 30.0008 
Y(ll)  a 28.01055 

Y C 12)  a 2.016 

Y C 13)  a 12.011 
YI14)  a 14.008 

C MOLECULAR  WEIGHT  OF  AIR 

MWA1R  a 28,9646 

C HYOROOEN/C ARSON  ATOM  RATIO 

MaO.1636 

R » 1.987165  * 778.28 

RA  a R/MVA1R 

R8  a (RA  . (IF  • R • M)  / 1 2,*Y  ( 12)  * (1.0  • •«>>>)  / (1.0  ♦ F) 

XMVGAS  a R/RB 
WTM  a XMWGAS 

C Check  FOR  ENTHALPY  KNOWN  case  to  calculate  TEMPERATURE 
IF (HH) 10002.10000,10001 

10001  TIa3000.*182.*ALOG10IPAI 

HI  »0 
SSaO 
CPP»0 

CALL  GASITI.PA.HI.SSiCPP.F.M.LP.ITEM) 
c EXHAUST  gas-check  for  oissocation  or  NON-OISSOCATION 
lF(HI-HH) 10002.10004,10003 
10004  TR  a T1 
RETURN 

10003  CALL  GAS(TR,PA,HM,SS,CPP.F,M,LP»ITEM> 

RETURN 

10002  ICONal.O 

TR  ■ 3000.0 

CALL  CHEMEO (PA, TR, F.HV, ICON, XKON.HH.SS.CPP.WTM.LP. ITEM) 

RETURN 

C CHECK  FOR  AIR  OR  EXHAUST  GAS 
10000  IF(F)5l, 100.50 

50  T I a 3000,0  • 182.  • ALOGIO(PA) 

IFITR.GT.TDGO  TO  59 
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FORTRAN  IV  6 LEVEL  21 


THERMO 


DATE  « T5224 


13/41/03 


0049 

0050 

0051 

0052 

0053 

0054 

0055 

0056 

0057 
0050 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0067 
0060 

0069 

0070 


CALL  OAS<TR.PA.HM,SS.CPP.F,M,LP.lTEM> 

RETURN 
TI»5000.0 
59  CONTINUE 
ICON  ■ I 
MR-0.0 
HG-0 

CALL  CHEMEO (PA. TR.F.HV. ICON, xKON,MH,SS.CPP.VTM,LP, ITEM) 

CALL  CM£MEQ<PA,TI,F,HV.!C0N,XK0N,MR,SR,CPR,*TR,LP,ITEM> 

Call  gasui.pa.hg.sg.cpg.f.m.lp.itemi 

OELM  - HH  - HR 

MM  - MG  • OELM 

OELCP  - CPP  - CPR 

CP P » CPG  * OELCP 

OELS  * SS  - SR 

SS  » SG  * OELS 

RETURN 

AIR  EOUATIONS 

100  CALL  GASITR.PA.MM.SS.CPP.F.M.LP.ITEM) 

WTM  ■ MWAIR 
DPTliON 

51  VRITE(LP.52)IT£M,F 

52  FORMAT ( IM  , I4,2X ,6MTM£RM0,2X, I 5M INCORRECT  INPUT. 3X.3MF  -.FI2.4) 
RETURN 

End 
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FORTRAN  IV  G LEVEL  21 


GAS 


OATE  « 7522A  13/41/03 


0001 

0002 

0003 

0004 

0005 

0006 
0007 

oooe 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 
0017 

ooie 

0019 

0020 
0021 

0022 

0023 

0024 


0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 

0039 

0040 

0041 

0042 

0043 

0044 

0045 

0046 

0047 

0048 

0049 

0050 

0051 


SUBROUTINE  GAS t TR ,PA,mB .5B.CP8.F ,H,LP. ITEM) 

REAL  M.MWAlR 
OATA  CTR2/0.0/ 

C0MM0N/CATMY/XN2,X02,XAR,XC02,XNE,MyAIR,Y(14> , CHASCO  < 1 0 ) » PAR » PNE 
OIMENSION  M<6),5<6>,CPI6>,Cl(4>,C2t4>,C3(4),Hl<72>,51<72>,CPl(72>, 
1X16) 

X(l)  ■ XN2 
X (2)  • X02 
X (3)  ■ XAR 
X(4)  ■ XC02 
X (5)  ■ XNE 
X (6)  ■ 0.0 
HIG*0.0 
MC«0.0 

lF(F.EQ.0.)HB-0.0 
IF (F.EO.O. ) GO  TO  74 
CALL  MFNODtX(l).F) 

IF (HB.GT. 1 00.0) G0T072 
IF(mB.EG.>1.0)GO  TO  72 

74  IFITR-4000,0) 75.300.999 

75  IFITR-300. )999«200.76 

76  IFITR-900. 01200.300.300 

C AIR  PROPERTIES  300R.-900R. 

200  CALL  MSCP (TR. 1.1. Mill) 

CALL  MSCP (TR. 1 .2.5(1)) 

CP(1)  ■ ((((-0.442I008E-12  *TR)  ♦ 1 .5962495E-09) *TR  - 

1 1 . 3S6283E-06)  • TR  » 4. 1982143E-04)  * TR  » 

2 6.9153738 
CALL  MSCP (TR ,2, 1 ,H (2) ) 

CALL  HSCP<TR.2.2.5(2>> 

CALL  MSCP (TR.2.3.CP (2>  1 
M C 31  ■ 4.96B15*TR 

5(3)  ■ 4.9681S  • ALOG ( TR/ 1.8)  ♦ 8.6765909 

CP (3)  « 4.96815 

CALL  MSCP(TR.4,1,M(4) ) 

CALL  MSCP (TR.4 .2.5 (4) ) 

CALL  MSCP (TR.4. 3. CP (4) ) 

M<5>  ■ 4.96815  • TR 

5(5)  ■ 4,96815  • ALOGITR/l .8)  • 6.6417222 

CP (5)  ■ 4,96815 

CALL  M5CP(TR.5.1.M(6) ) 

CALL  M5CPITR .5. 2, 5(6) ) 

CALL  MSCP (TR .5. 3. CP (6) ) 

DO  6 1*1.6 

IF  < X ( I ) .LE.O.IGO  TO  5 

Sill  ■ 5(1)  >1.987165  • ALOG ( X ( I ) . • PA) 

GO  TO  6 

5 SU)  ■ 0. 

6 CONTINUE 
GO  TO  1000 

C AIR  AND  REVISED  PROPERTIES  900R-4000R 
72  MC»0 
MC»MB 
TR  • 2Sr/0. 

IF(MC) '04,104,300 
300  CALL  CCEFF (7, Ml) 


59 


AEDC-TR-76-15 


! 


FORTRAN  IV  0 LEVEL  21 


GAS 


DATE  « 7522* 


13/*l/03 


0052 

CALL  COEFF (8.S1 ) 

0053 

CALL  COEFF (9.CP1 ) 

005* 

IF(TR  .L£.  *000.0)  K A**8 

OOSS 

IFITR  ,L£.  2700.0)  KA»2* 

0056 

IF  ( TR.LE « 1800. ) KA*0 

0057 

DO  *00  K*1.6 

0058 

H(K>  ■ 0.0 

0059 

S(K)  ■ 0.0 

0060 

CP ( K > * 0.0 

0061 

KK  * * •(«•!)  «KA 

0062 

DO  500  J«I.* 

0063 

KKK  * KK.J 

006* 

TX  ■ T9  ••(J-I) 

0065 

CKJ)  ■ HI  (KKK ) • TX 

0066 

C2 ( J)  ■ SI (KKK)  • TX 

0067 

soo 

C3(J)  • CPI (KKK ) *TX 

0068 

DO  30  J*I.* 

0069 

H(K)  ■ H(K)  * Cl (J) 

0070 

S(K)  ■ S(K)  * C2IJ) 

0071 

30 

CP (K)  ■ CP(K)  « C3< J) 

0072 

IF(K*I ) 320.315.320 

0073 

315 

IFITR  ,LE.  2700.0. AND.  TR  .GT.  1800.0ISIK) 

■ S(K) 

I • TR  ••  * 

007* 

320 

CONTINUE 

0075 

*00 

CONTINUE 

0076 

DO  IS  1*1.6 

0077 

IF  (X(I).LE.O.IGO  TO  1* 

0078 

S(I)  • Sill*  I. 987165  • ALOGCXII)  • PA) 

0079 

GO  TO  IS 

0080 

I* 

S ( I ) * 0. 

0081 

15 

CONTINUE 

0082 

1000 

MA  * XN2*  MID  «X02*H ( 2 ) .XAR  *H(3)  » XC02 

•H(*) 

0083 

HAIR  * HA/MVAlR 

C 

CALCULATION  OF  PROPERTIES-AE0C-TR-7S- 

-SECT 

008* 

AH  . (H (6)  -.5  * H (2) » / Y 1 12) 

0085 

BH  « (H(A)  - H (2) 1 / Y ( 131 

0086 

2 « l./Il.»F) 

0087 

HB  • 2 * (HAIR  . F *(AH  • H .BH)  /(l,  «H) 1 

0088 

SA  * XN2  ‘SID  ♦X02«S(2)*XAR*S(3I»XC02*S(*)»XN£  • 

0089 

SAIR  * SA/MUA1R 

0090 

AE  ■ (S (6)  - .5  • 5(2)1  / Yf  121 

0091 

BE  • (SI*)  - 5(2)1/  Y ( 13) 

0092 

SB  ■ 2 • (SAIR  ♦ F* (AE*H  ♦ BE)  / (1.  ♦ HI) 

0093 

CPA  ■ XN2  «CP<1)  ♦X02*CP(2)*XAR*CP(3) •XC02*CP(*) 

009* 

CPAIR  * CPA/MVA1R 

0095 

AC  « (CP(6)  » .C  • CP ( 2 ) ) /Y(12) 

0096 

BC  *(CP(*)  - CP (21 ) /V ( 131 

0097 

CPB  • (CPAIR  ♦ F • ( AC*H  •BC)/(H*1.))  / (1. 

♦F) 

0098 

71 

IF(HC  .GT.  i.E-6)  GO  T073 

0099 

10* 

IF (F) 105.150.105 

0100 

73 

CELT  * (HC-HBI/CPB 

0101 

IF  (ABS(DELT) .LE.0.001*TR)GOT0105 

0102 

CTR2  • CTR2M.0 

0103 

IF (CTR2.GT .25.0) GO  TO  10S 

010* 

lF(OELT)  70.105.70 

0105 

70 

TR*TR»0ELT 

H (51 


2.2.1 


• S (SI 


XNE  «CP(51 


I 
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0106 

0107 

0108 

0109 

0110 
0111 


0112 

0113 

0116 

0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 


FORTRAN 

0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 


IV  G LEVEL  21  GAS  OA7E  • 75224  13/41/03 

GO  TO  300 
69  TRaTR-OELT 

WRITE(6.68)TI.PA,H1G.MC,HB 
68  FORMAT! 12H  GAS  ■•10EI2.6) 

GO  TO  300 
105  RETURN 

C CALCULATION  OF  AIR  PROPERTIES  WITH  INTEBMOLECULAR  CORRECTION 
C CALCULATION  OF  PROPERT 1ES-AE0C-TR-75-  ‘SECT. 2.1. 4 

150  CPB  ■ CPMIXP(PA.TR.CPAIR.LP.ITEM) 

H8  ■ HMlXPIPA.TR.HAlR.LPi  ITEM) 

SB  ■ SMI XP (PA (TR. SA1R.LP. ITEM) 

return 

999  WRITE (LP.250) ITEM 

250  FORMAT (IM  . 14.2X.6HTHERM0.2X.25H  TEMPERATURE  OUT  OF  RANGE) 

H8  ■ 0.0 
SB  ■ 0.0 
CPB  ■ 0.0 
RETURN 
END 


IV  G LEVEL  21  MFNOO  DATE  ■ 75224  13/41/03 

SUBROUTINE  MFNOO(X.F) 

C CALCULATION  of  EXHAUST  GAS  C0MP0S1TI0N-T.LE.TI, 

REAL  MWA1R,MN2.M02,MAR,MC02.MNE.MH20.MWFL 

COMM ON/CA THY/ XN2. X02.XAR. XC02. XNE.MWA1R.V ( 14) ,2(10) .PAR.PNE 

OIMENSION  x (6) 

XN  - 1.0 
XM« 1 .95 

MwFL  ■ 12.011  • XN  . 1.00797  • XM 
FAM  a F • MWA1R/MWFL 
MN2  » XN2/FAM 

M02  * X02/FAM  - XN  - XM/4. 

MAR  ■ XAR/FAM 
MC02  a XC02/FAM  ♦ XN 
MNE  a XNE/FAM 
MH20  a XM/2. 

TM  a MN2  »M02*MAR  »MC02  • MNE  ♦ MM20 

X < 1 ) a MN2/TM 

X (2)  a M02/TM 

X (3)  a MAR/TM 

X (4)  a MC02/TM 

X (5)  a MNE/TM 

X (6)  a MH20/TM 

RETURN 

ENO 
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0001 

0002 

0003 


000* 

0005 

0006 

0007 

0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 


? 


f, 

i 


G LEVEL  21  CMENEO  OATC  • 75224  13/41/03 

SUBROUTINE  ChEMEQ(PA,T.F,MV,1C0N,XK0N,M,$,CP.*TM,lP,1TEM) 

REAL  M 

C0MM0N/CATHY/XN2, X02  «XAR<  XC02 . XNE, MVA 1R»  Y 1 1 4) • X ( 1 0) ,PAR,PNE 

c this  subroutine  initializes  the  rmooes  ano  osgerby  eouations 

C TO  CALCULATE  THE  PROPERTIES  OF  A OISSOCUTEO  EXHAUST  GAS 
XN  « 1.0 
XH  • 1.95 

XX  • 0.209476/ ( XN  « 0.25  • XH) 

C STOICHIOMETRIC  RATIO  * 0.06795  fOR  CNH1.95N 
PHl  a F/  0.06795 
C a 2.0  • X02  * 2.0  • XC02 
0 a 2.0  • XN2 
E a XAR 

GaxNE 

A a XC02  ♦ PHI  • XX 
8 a PHI  • XX  • XH 
1F(XK0N) 10.20.10 

10  Call  fuElhia.b.phi.xx.xm.hv.h) 

GO  TO  66 
20  CONTINUE 

IF(h.GT.O..ANO.XKON.EO.O.)GO  TO  30 

CALL  PROP(A.B.C.O.E«G.T.PA.H.S.CP.*TM,x.PAR.PNE.TM,NT.1CON.LP. 
S1TEM) 

RETURN 

30  TIa3000. * 182.*AlOG1 0 (PA ) 

Ma ( XM«Y ( 7)) /(XN»Y( 1311 

CALL  r.AS(Tl.PA,HIG,SS.CPP.F,M,LP.lTEM) 

DELH  a H • HlG 
HIPaO.O 

CALL  PROPU.B.C.O.E.G.TI  .PA.HIP.S.CP.HTM.x.PAR.PNE.TH.NT.  ICON, 
1LP.1TEM) 

H a HIP  « OELH 
66  ICONao.O 
TaO.O 

CALL  PROP(A.B,C.O.E.G,T.PA.H,S,CP,HTM,X,PAR,PNE.TM,NT,ICON, 
1LP.ITEM) 

RETURN 

ENO 
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0001  SUBROUTINE  PROP (CAR.HYD. 0* • AN2.AR.ANE. TT.PR.H, S.CP.WTM.P.PAR.PNE t 
•TIHE.NT.IT.LP.ITEM) 

0002  C0MM0N/CATHY/XNT,X02.XAR.XC02.XNE.MVAIR,Y(I41 

C SPECIES  ARE  I-H, 2-0, 3-0H.A-N0.5-CD.6-H2, 7-02, 8-H2O.9-CO2.10-N2. 

C I1-N.12-C.1I-N.I2-C 

C SUBROUTINE  FOR  DETERMINING  EQUILIBRIUM  PARTIAL' PRESSURES 

C T IS  TEMP.  IN  OEG  R ANO  CAN  BE  EITHER  INPUT  OR  OUTPUT 

C P IS  PRESS  IN  ATMOSPHERES  ANO  IS  AN  INPUT 

C H IS  THE  ENTHALPY  IN  8TU  PER  LBM  ANO  IS  EITHER  AN  INPUT  OR  OUTPUT 

C S IT  THE  ENTROPY  IN  BTU  PER  LBM-OEG  R ANO  IS  AN  OUTPUT 

C CP  IS  THE  SPECIFIC  HEAT  IN  BTU  PER  LBM-OEG  R ANO  IS  AN  OUTPUT 

C WTM  is  THE  MOLECULAR  *T  OF  THE  MIXTURE  ANO  IS  AN  OUTPUT 

C NT  IS  THE  NO.  OF  ITERATIONS  TO  FINO  PARTIAL  PRESS  ANO  IS  AN  OUTPUT 

C IT  IS  A CONTROL  CONSTANT  FOR  THE  SUBROUTINE 

C ITbO  h »N0  P ARE  KNOWN  *N0  T.CP.S.  ANO  WTM  ARE  FOUNO 

C *T«1  T A P ARE  KNOWN  • H.S.CP.WTM.ANO  PARTIAL  PRESSURES  FOUNO 


0003 

DIMENSION  01  (3)  .D2C 

•(3). 03(3! -C (61 

0004 

DIMENSION  DAR (31 .0NEI3) 

0005 

OIMENSION  F (121  .A1  (3)  •/ 

C 

CHANGE  units  of  inputs  ' 

0006 

T ■ <5./9.1»TT 

0007 

H0«H 

0008 

C 

ICT«0 

0009 

OELTLaO. 

0010 

00  300  1«1.I0 

0011 

300  P(I1«0. 

0012 

DELL*. 0001 

0013 

FUEL*2.«CAR*.S»HY0 

c 

RATIO  OF  ELEMENTS 

0014 

RAN  ■ AR/AN2 

0015 

RCNm  C AR/AN? 

0016 

RHNbHY0/AN2 

0017 

RONbOX/ANZ 

0018 

/» 

RNEN* ANE/ AN2 

0019 

V. 

0£LI»0. 

0020 

0EL2»0. 

0021 

DEL3«0. 

0022 

IF (UX.GE .FUEL  1 GO  TO  l 

0023 

WRITE (LP .971  ITEM 

0024 

FUEL  b DX 

0025 

1 CONTINUE 

c 

initial  guess  of  o2,H2o 

0026 

XX  S ,5/t (.S*HY0*DX«AN2 

0027 

P (7) bPR»xX» (OX-FUEL) 

0028 

P(9)b?,»CAR*XX»PR 

0029 

P(8)bHYO*XX»PR 

0030 

P ( 1 0 1 * AN2»XX»PR 

0031 

PaRb2,#P(101«RAN 

0032 

PNE«2,  • P(10)  • RNEN 

0033 

L«  7 

c 

ITERATION  LOOP  TO  FIND  1 

0034 

IF  (P(L1.LE..011  P(L>».' 

0035 

00  4 Kb] .SO 

..  •»  a sO££,  Satofellinifaiifc 
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0036 

C 

KOUNTbK 

CHECK  all  initial  partial  pressures  FOR  NEGATIVE  values 

0037 

00  41  jBl.10 

0030 

C 

41 

IP  (P(J).LT.O.)  P ( J) ■ 1 «E*5 

OON* T RECALCULATE  FREE  ENERGV  RELATION  TOR  T KNOWN— 1T«1 

0039 

IP  (K.GE.2.AND.1T.EQ.1)  GO  TO  10 

0040 

IF<1T.EO.O.ANO.K.EQ.1>Tb16TO. 

0041 

CALL  TH£RM(P,T.M,5,CP,F.PR.PAR,PNE,WTM,0) 

0042 

E(l)  » 2.71828182flS»«(.S»PIS)-.25»P!7)-P(l)) 

0043 

E(2)  ■ 2.71820182flS»«(.S»P(7)-P(2)) 

0044 

E 13)  ■ 2« 71828182*5  ••  I.S  • P (8) ».25»P 1 7) -F ( 3) ) 

004S 

£(4>  • 2.7182818285  ••  (.5  •FI10)  • .5  • F(7)  -F(4)) 

0046 

£(S)  ■ 2.7182810285  ••  IP  <9)  - ,5  • P 1 7)  - P(5)> 

0047 

C 

E (6)  » 2.7102018285  ••  IP  IS)  - .5  • r 17)  - F(6)l 
CALCULATE  PARTIAL  PRESSURES  OP  THE  ASSUMED  KNOWN  SPECIES 

0040 

10 

Pills  SORT (P (0) /SORT (PIT) ) )•£( 1 ) 

0049 

P (2)  ■ SORT  IP  1 7) ) •£ 12) 

0050 

P<3)«SQRT<P<0) ‘SORT (PIT) I l»E(3) 

00S1 

P <4IBS0RT<P< 10) »P (7 ) )•£ <41 

0052 

XXal ./ ( 1 ,♦£ 151/SQRT IP (7) ) ) 

0053 

P<9)b  RCN*  < 2«*P  < 10) *P<4) )*XX 

00S4 

P < 5) ■ P< 91 /SORT  IP (7) ) bE(5) 

0055 

P<6|b  PI81/S0RT (PIT) ) *E<6) 

00S6 

PAR  a (2,  s p(lo>  • P ( 4) ) BRAN 

0057 

PNEb  (2.  b P(10)  ♦ P (4) 1 • RNEN 

w 

c 

DERIVATION  OP  RELATIONSHIP  OP  KNOWN  SPECIES  WITH  RESPECT  TO 

0050 

00  11  1b1.3 

i.059 

Olltlao. 

0060 

02 111 a0 • 

0061 

031 1)  «0. 

0062 

04(11 «0. 

0063 

OSIllBO. 

0064 

06<  n«o. 

00  65 

0711)80. 

0066 

00(1)80. 

0067 

09<1)bO. 

0060 

11 

oio<!ibo. 

0069 

01111b  -P(l)/p<7l»,25 

0070 

01121b  P< 1 ) /P (8) *.S 

0071 

02<1)b  PI2)/P(7|b,5 

0072 

03(1)b  P(3I /PI  71 b.2S 

0073 

03 (21 ■ P(3)/P(8)»,5 

0074 

04(1)b  P(4)/P(T|b,5 

0075 

04 ( 31 »P( 41 /P ( 10) *.5 

0076 

IF  (CAR. EO. 0.01  GO  TO  441 

0077 

09(1)8  RCN»(04(1Ib(2.»P(10)»P(4))»XX».S*E(5)/  PIT) 8»1 .5)»XX 

0070 

09(3)b  RCNBXX*(2.’04(3)) 

0079 

05(1 )■  P(5)/P(9I»09(1)-P(5)/P(7)*,5 

0000 

05(31b  P(5)/P(9I»09(3) 

0001 

441 

CONTINUE 

0002 

06(1)b  -P(6I/P(T)B,5 

0003 

06(2)b  P(6I/P(0) 

0004 

07<llal.O 

OOOS 

00(2)si,0 

0006 

O10(3)al,0 
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0087 

0088 

0089 

0090 

0091 

0092 


0093 

0094 

0095 

0096 

0097 
009P 

0099 

0100 
0101 
0102 
0103 


0104 

0105 

0106 

0107 

0108 


0109 

0110 
0111 
0112 

0113 

0114 


0115 

0116 

0117 

0118 

0119 

0120 
0121 
0122 

0123 

0124 

0125 

0126 
0127 


OAR I 1 ) *RAN*04 ( 1 1 
OAR  ( 2)  * 0. 

OAR  t 31 *RAN* (04(3) *2*1 
ONE  1 1 ) • RNEN  « 04(11 
ONE (2)  ■ 0.0 

ONE ( 31  * RNEN  > (04(3)  * 2.) 

C 

C DERIVATION  OF  MATRIX  coefficients 
0012  1*1.3 

AI(I)«  Dld)*03d)*2.*D6<I)*2.*D8(I>-»MN*(2.»DlO(I)*D4ll)) 

A2(I>«  02  (11*03(1)*  04  d)*05d)*2.*07(I)*  08(1) *2. *09 ( 1 1 “RON* (2. 

**0I0(I)*04(I)) 

A 3 ( I ) * 01 (I) *02(1) *03(11*04(1) *05(1) *06(1) *07(11*08(1) *09(1) *010 1 
*I)*DAR(I)  ♦ ONE (I) 

12  CONTINUE 

SP1  « P(I)*P(3)»2.#P(6) *2.»P  (8)  -RHN* (2.*P  1 10) »P(4) ) 

SP2  ■ P( 21 *P (3) *P (4) *P (5) *2.»P(7)*P(8)*2.»PC9)-RON»(2.»P(10)*P(4 

•I  I 

SP3  *0. 

00  13  1*1.10 
11  SP3  «SP3  *P(I) 

SP3«SP3-PS«PAR  • PNE 

C 

C SOLUTION  OF  MATRIX 

OEN*  A1 n )« (A2(2)*A3(3)-A3(2)*A2(3) )-A2(l)*(Al (2) *A3 ( 3) -A3 ( 2) • 

•A  1(3) I *A3( I)»(A1 (2)«A2(3)-A2(2)*A1 (3) ) 

OEN* 1. /OEN 

XNI  *SPl*(A3(2)*A2(3)-A2(2)*A3(3)) *SP2* (A1 (2) *A3t3) -41 (31 *4342) ) * 
•SP3*(A1 (3)*A2(2)-AI (2) *A2 (3) ) 

XN2  «SP1*(A2(1)*A3(3)-43(1)*A2(3) !»SP>»(A3(l)»Al(3).*l(l)4*3t3>)* 
•SP3*(AI (1)*A2 ( 3 ) — A 2 ( 1 ) • A I (3) ) 

XN3  *SPt»(A3(l)*A2(2)-A2(D*A3(2)  )*SP2*(AHI)»A3(2)-A3(1)*A1(2)  )♦ 
•SP3* ( A2( 1 ) »A1 (21-A1 (ll»A2(2) ) 

C 

C CALCULATION  of  PARTIAL  PRESS  CORRECTION  AND  NEW  GUESSED  PART  P»ES5 
OEL 1 * XNI  *OEN 
0EL2*  XN2  »OEN 
0EL3*XN3*0FN 
P(L)*P(L)  *0v.-' 

PIS) «P (81 *0EL2 
P ( 10>  *P( 10 ) *0EL3 
C 

C FIND  ENTHALPT  FOR  ITERATION  TO  FINO  temp. 

7 CONTINUE 
IF (IT) 18.18.6 

18  CALL  THERMIP.T.M  .S.CP.F.PR.PAR.PnE.WTM.I  ) 

80  CONTINUE 

IF  (1CT.GT.0)  GO  TO  85 
OELT  • ((M0-HI/CP1/1.8 
IF  (OFlT*OELTL)  81.81. 82 
82  T»T*0FLT 
GO  TO  85 

81  T.T.OELT 
ICT-2 

85  1CT«ICT-1 
0ELTL«0ELT 
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0128 

0129 


0130 

0131 

0132 

0133 
013* 

0135 

0136 


0137 


0138 

0139 
01*0 
01*1 

01*2 


1F<T.LE.300.)T*300.0 
IF  (T.GE.5000.)  T*5000. 

C 

C ITERATION  CONTROL— TEMP. VI THIN  I 0E6  OR  PRESS  VlTHlN  DELL 
IF ( ABS (HO-H)  -.1*  CP) 6.6.* 

6 IF ( ABStOEL  1 ) -<O.OI«OELL> > 22.22.* 

22  IF  (ABS(0EL2)  - DELL)  23.23.* 

23  IF (ABSIDEL3)  - DELL)  5.S.* 

* CONTINUE 

5 CONTINUE 
PHOLOaPR 

C VR1TEILP.3) ( P ( I > , 1-1*10) 

C 3 FORMAT 1 1 OFT.*) 

CALL  THERMlP.T.M.S.CP.F.PR.PAR.PNE.VTM.l) 

C WRITE (LP.3) <P( I) • 1*1,10) 

C 3 FORMAT ( 10F7.*) 

TT*T*1 ,8 
NTaKOUNT 
RETURN 

97  FORMAT ( IN  ,I*,2*.6MPR0P  .2X.36H  THIS  CASE  HAS  TOO  LITTLE  02  TO  RU 

•N  ) 

ENO 


66 


Afli1  ■Tl’MiiWn 


AEDC -TR-76-15 


FORTRAN  IV  G LEVEL  21 


OATE  ■ 7522* 


0001 

SU8R0UTINE  THERMIBLP.T.M.S.CP.F.PR.PAR.PNE.WTM.i 

0002 

COMMON/CATHY/XN2. X02. XAR .XC02.XNE .MVA1R. Y ( 1 4) 

0003 

DIMENSION  ALP (12) .CX (72) «CY(72) .CTTI144) .HXI12) 

1C 1 (4) .C2 (4) <CT (72) .F ( 10) «CF (4) .C3(4) 

0004 

DIMENSION  8LP (10) 

0005 

EQUIVALENCE (CTT ( 1 ) .CT ( 1 ) ) . (CTT ( 73) .CX ( 1 ) ) 

0006 

00  50  1>1.6 

0007 

50 

ALP ( 1 ) ■ 8LP ( 1 ) 

0006 

ALP (7)  ■ BLP(IO) 

0009 

ALP (8)  ■ 8LPI7) 

0010 

ALP (9)  ■ PAR 

0011 

AlP (10)  ■ BLP (9) 

0012 

ALP (11)  ■ PNE 

0013 

ALP (12)  a BLP (8) 

0014 

00  51  lal . 12 

0015 

51 

ALP(l)  a ALP ( I ) / PR 

0016 

CPaO.O 

0017 

H a 0.0 

0016 

5 a 0.0 

0019 

1FIKT) 40 0.100. 200 

C 

CALCULATE  FREE  ENERGY 

0020 

100 

CALL  COEFF(l.CT) 

0021 

T a T • 1.6 

0022 

00  201  Kal.lO 

0023 

F(K)  a o.o 

0024 

KFa4» (K«l ) 

0025 

00  300  Jal.4 

0026 

KFF  a KF  • J 

0027 

TX  a T aa(J-l) 

0026 

300 

CF(J)aCT(KFF)*TX 

0029 

00  301  Jal,4 

0030 

301 

F(K)  a F (K)  * CF(J) 

0031 

201 

CONTINUE 

0032 

F10  ■ F ( 10 1 

0033 

F ( 10)  a FIT) 

0034 

F (7)  a F (8) 

0035 

F ( ft)  a F10 

0036 

00  81  I>1. 10 

0037 

61 

Fill  a (F  ( 1 ) /T)» (1 .6/1 .967165) 

0036 

TaT/1.8 

0039 

RETURN 

C 

CALCULATE  HICHEMICAL  ♦ SENSIBLE  ENTHALPY) .5. CP 

0040 

200 

CALL  COEFF (2.CT) 

0041 

CALL  COEFF ( 3.CX) 

0042 

CALL  COEFF (4, CT) 

0043 

T a T»1 .8 

0044 

00  2 Kal . 12 

0045 

HX(K)  a 0 

0046 

SX(K)  a 0 

0047 

CPX(K)  a 0 

0046 

KA  a (| 

0049 

IF (K>7) 20.21 .21 

0050 

21 

IFIT-2700. 0)20.20.27 

0051 

27 

KA  a 24 

0052 

20 

KK  a 4* (K«l ) • KA 

0053 

00  3 J«  1 .4 

AEDC  TR-76  15 


FORTRAN 

0054 

0055 

0056 

0057 

0058 

0059 

0060 
0061 
0062 

0063 

0064 

0065 

0066 
0067 
0069 

0069 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 


0079 

0080 
0081 
0082 

0083 

0084 

0085 


FORTRAN 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0011 


IV  6 LEVEL  21  THERM  OATE  ■ 75224  13/41/03 

KKK  ■ KK«J 
TX  a T ••Oil 
CKJI  ■ CT(KKK)  • TX 
C2(J)  ■ CXI KKK I • TX 
3 C3( J)  ■ CY (KKK I • TX 
00  30  J*1 *4 

HX (K I a HXIKI  ♦ C1CJI 
Sx  (K ) a SX(K)  ♦ C2(  J) 

30  CPX (K | ■ CPX(K)  ♦ C3( J) 

1F(ALP(KI.L£. 0.160  TO  11 

SX (K ) ■ SX(K)  - 1.987165  • ALO0(ALP(K)  • PR) 

60  TO  12 

11  SX (K) »0. 

12  CONTINUE 
lF(K-6)25. 23.25 

23  CPX (K I a CPX (K)  ♦ 0.439912E-I3  *T  «*4.0 
25  lF(K-7)35. 29,35 
29  IF(T>2700, 0)31,31, 35 

31  SXCK)  a SX (K ) ♦ 0.6848142E-12  • T ••  4.0 
35  CONTINUE 

H a H « ALP(K)  • HX  (K) 

S a S . ALPIK)  • SX(K) 

CP  a CP  ♦ ALP(K)  • CPX(K) 

2 CONTINUE 

VTM  a ALP(l)  V(7)  ♦ ALP(2>  • Y<8)  * ALP (3)  • Y(9)  • ALPC4)  • 

1 V(1S)  • ALP15)  • Y ( 1 1 1 ♦ ALP (6)  • YI12)  ♦ ALP (7)  * Ytl)  ♦ 

1 ALP (8)  • Y(2)  • ALP (9)  • Y(3)  • ALP(10)*Y(4)  ♦ ALP(ll)  * Y(5) 

1 ♦ ALP (12)  • Y<6)  

H a H / VTM 
S a S/MTM 
CP  a CP/WTM 
T a T / 1.8 
400  CONTINUE 
RETURN 
END 


IV  6 LEVEL  21  FUELH  OATE  a 75224  13/41/03 

SUBROUTINE  FUELH(A.B.PHI,XX,XM,HV,M) 

C CALCULATES  ENTHALPY  OF  FUEL  FROM  HEATING  VALUE. 

C0MM0N/CATHY/XN2,X02.XAR,XC02,XNE,MVA1R,Y(14) , DUMMY! 12) 

C FUEL  ANO  AIR  ARE  AT  536  DEGREES  R 
TFUEL  a 536.7 

C EO.  FOR  CPFUEL  FROM  NASA-TN-3276 

CPFuEL  a 0.0005* (TFUEL-459. 671 »0. 44 
HAIR  ■ >1.8680522 

PCTF  a ( Y ( 13)  • A • Y(7I  • 8)/(Y(13>  • A ♦ Y(7)  • B ♦ Y(8)  • 2.0  • 

1 X02  ♦ Y (14)  • 2,0  • XN2  ♦ Y<3)  • XAR  ♦ Y(5)  • XNE) 

VTMF  a y( 13)  • A ♦ Y(7)  • B 

HF  a MV  - (94054,0  • (A--.C02)  ♦ 57797.9  ♦ (B/2.0) ) •<  1 .8/VTMF) 

• ♦ CPFUEL  • (TFUEL-536.V) 

H a PCTF  • MF  ♦ (l.o  - PCTF)  • HAIR 

RETURN 

ENO 


68 


naaaim 


J 
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HMIXP  OATE  ■ T522A  13/41/03 


0001 

0002 

0003 


0004 

0005 

( 0006 

0007 

! 0008 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 


0025 

0026 
0027 


FUNCTION  HNIXP(PA«TR«HIO,LP, ITEM) 

REAL  MWAIR 

CONNON/CATMY/ANF (5) .MWAIR, Y ( 14) 

C CALCULATE  ENTHALPY  WITH  PRESSURE 

C CORRECTION 

CH  ■ 82.056 

C GAS  CONSTANT(ATM-CM**3)/(GM-MOLE-(OEG)K> 

T * TR/I.8 
IF(T-166. 13.4,4 

3 T*I66.0 

4 OHM  * 0.0 

OELT  ■ (CH  • T / PA)  / 10000.0 
CTR  ■ 0.0 
00  20  1*1 .5 
ei  * RBB(I.I.TR) 

Cl  • CCC(I.I.TR) 

BPI  ■ 866(2,1, TR) 

CPI  * CCCI2.I.TR) 

C CALCULATION  OF  SPECIFIC  VOLUME! !> 

VIG  ■ CH  • T / PA  * BI  ♦ Cl 

14  VI  * CH  • T / PA  • (1,0  • BI  / VIG  ♦ CI/IVlG  • VIG) ) 

IF ( A6S( Vl»VlG) »OELT> 18,18,15 

15  CTR  ■ CTR  ♦ 1,0 
VIG  ■ Vl 

IFICTR-25.0) 14,16,16 

16  VI  * CM  • T / PA 
18  CONTINUE 

OH  a 1.987165  / Yd)  • 

I TR»  ( 1.0  / VI  • (BI  - BPI)  » I. 0/VI  • I. 0/VI  a(CI  - 0.5« 

2CPI ) ) 

20  OHM  a OH  • AMF(I)  « OHM 
HMIXP  a MIO  • OHM 
RETURN 


FORTRAN  IV  G LEVEL  21 


SM1XP  DATE  ■ 7522*  13/41/03 


i 

t 


0001 

0002 

0003 


000* 

0005 

0006 
0007 
0006 

0009 

0010 
0011 
0012 
0013 
001* 

0015 

0016 
0017 
0016 

0019 

0020 
0021 
0022 
0023 
002* 


0025 

0026 
0027 
0026 


FUNCTION  SMlXPIPA.TR.SIO.LP.ITEK) 

REAL  MVAIR 

COMMON/C ATHY/AMF (5)  .MVAIR.Vd*) 

C CALCULATE  ENTROPY  WITH  PRESSURE 

C CORRECTION 

CS  ■ 62.056 

C GAS  CONSTANT  ( ATH«CM«»3) / (GM-MOLE- < OEGI XI 

T ■ TR/1.6 
IF(T-166.)3.*»* 

3 T«I66.0 
* DSN  ■ 0.0 

OELT  a (CS  • T/PAI  / 10000.0 
00  20  I a 1.5 
B1  a BBB(l.I.TR) 

Cl  a CCC(l.l.TR) 

BPI  a SBB (2*1  * TR 1 
CPI  ■ CCCI2.I.TR) 

CTR  a 0.0 

C CALCULATION  OF  SPECIFIC  VOLUME(l) 

VIG  a CS  • T / PA  ♦ BI  » Cl 

I*  VI  ■ CS  • T / PA  • (1.0  *81/  VIG  ♦ Cl  / (VIG  • VlG)  1 

lFIABS(Vl-VlG)-DELT) 16*16.15 

15  CTR  ■ CTR  * 1.0 
VIG  a vl 

IF (CTR-2S.0) 1**16.16 

16  VI  a CS  * T / PA 
18  CONTINUE 

OS  a-1. 987165  / Y < 1 1 • 

1 < BPI/Vl  ♦ (1.0/  VI)  • (I.O/VI)  • ( 

2(BI  * Bll/  2.0  - CI/2.0  * CPI/2. Oil 
20  OSM  a OSH  * OS  • AHF ( I ) 

SMI XP  a SIO  ♦ OSM 

RETURN 

ENO 


AEDC-TR-76-15 


FORTRAN  IV  G LEVEL  21 


CPMIXP  DATE  ■ 7522*  IJ/M/03 


i] 

1 


0001 

0002 

0003 

0004 

0005 

0006 

0007 

0008 

0009 

0010 
0911 
0012 
0013 
001* 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 
0023 
002* 

0025 

0026 


0027 

0028 

0029 

0030 


FUNCTION  CPMIXP (PA.TR.CPID.LP. ITEM) 

C CALCULATION  OF  SPECIFIC  HEAT  WITH  PRESSURE  CORRECTION. 

COMHON/CATHY/AMFI5) «MVAIR.Y(t*> 

REAL  MmAIR 
Cl  ■ 82.056 

C GAS  CONSTANT  Cl (ATM-CM««3)/GM-MOLE  -<05G)K 

T « TR/I.8 
IF(T-|66.)3.*»* 

3 T«I66,0 
* OCP  *0.0 

CELT  « (Cl  • T / PA  ) /I0000.0 

00  20  I«I.S 

8PI  ■ B8B (2.1 «TRI 
CPI  • CCCI2.1.TR) 

81  ■ 888(1. I. TR) 

Cl  • CCC(I.l.TR) 

CPPI  ■ CCCI3.I.TR) 

BPPI  ■ 888(3. I. TR) 

CTR  ■ 0.0 

C CALCULATION  OF  SPECIFIC  VOLUME (I) 

VIG  ■ Cl  • T / PA  ♦ 81  » Cl 

I*  VI  ■ Cl  • T / PA  • (1.0  • 81  / VIG  ♦ CI/IV1G  • VIG) > 

IF  ( ABS( VI«VIG) -CELT) 18.18.15 

15  CTR  ■ CTR  « 1.0 
VlG  « VI 

IFICTR-25. 0)1*. 16.16 

16  VI  ■ Cl  * T / PA 
18  CONTINUE 

OCPI  « (-1. 987165)/  Yd) 

1 . * ( BPPI/  VI  - I. 0/VI  • I.O/Vl  * ((81  - BPI)  • 

2 (81  - BPI)  • Cl  • CPI  - 0.5  • CPPI) > 

20  OCP  « OCP  • OCPI  * AMFU) 

CPMIXP  • CPIO  ♦ OCP 

RETURN 

ENO 


f 

i 


i 
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r OUTRAN 
0001 

0002 

0003 

000* 

0005 

0006 


0007 


0008 


IV  G LEVEL  21  HSCP  DATE  ■ 7522*  13/41/03 

SUBROUTINE  HSCP(TX«K,KK.M) 

C CURVE  rn  ROUTINE  FOR  CONSTITUENT  PRDPERT1ES-300R-900R 

C - AEOC  -TR  -75-  -REF, 20 

01MENS10N  Cl <5,52,*) , C2(S.1*.*>.  C3<3,11.*) 

DIMENSION  XI  (95) , X2 (9S ) < X3 (95) «X* (95) .XS(95> ,X6(9S> .X7(95)  »XS(9S>, 
1X9(95),X10(95) ,X11 (90) «X12 (95) .X13 (95) ,X1* (90) «X15(95) ,X16(37) 
EQUIVALENCE  (X1(1).C1(1)),(X2(1),C1(96)).(X3(1),C1(191)). 
l(X4(l).Cl(286)>.(X5(l).Cl(381>>.<X6(l>,Cl<*76)),(X7<li.Cl(S71>>, 

2 <X8( 1 ) ,C1 (666) ),(X9(1),C1(761)).(X)0(1)«C1(B56)). 

3 (XI  1(1)  .Cl  (951) ) . ( X 1 2 ( 1 ) ,C2(1> ) ,<X13(1) »C2<96> ) • (XI*  < 1 ) ,C2(19D  ). 
*(X15(1) ,C3(1) ) • (X16 ( 1 ) ,C3 (96) ) 


REAL  INC 
Data  x i/ 


0.69300E 

00, 

0.69400E 

00, 

0.0 

• 

0.12518E 

04, 

0. 78600 E 

00, 

0.76300E 

00, 

0.76300E 

00, 

0.90000E 

01  9 

0.13777E 

04, 

0.86600E 

00, 

0.83300E 

00, 

0.83300E 

00, 

0.50600E 

03' 

0.15040E 

04, 

0.94600E 

00, 

0.90200E 

00, 

0.90300E 

00, 

0.10030E 

ot, 

0.16312E 

0*. 

0.10250E 

01. 

0.97200E 

00, 

0.97200E 

00. 

0.15000E 

04' 

0, 17594E 

04, 

0.11050E 

01. 

0. 10*1 OE 

01, 

0.10420E 

01, 

0.19960E 

04* 

0.1 888 8E 

04, 

0.11840E 

01, 

0.11110E 

01, 

O.llllOE 

01, 

0.24930E 

04' 

0.20197E 

04, 

0.12640E 

01. 

0.11800E 

01, 

0. 1 181 OE 

01, 

0.29900E 

04' 

0.21522E 

04, 

0.13440E 

01, 

0.12500E 

01, 

0.12500E 

01, 

0.34870E 

04' 

0.22864E 

04. 

0.14230E 

01. 

0.13190E 

01, 

0.13200E 

01 , 

0.39840E 

04' 

0.24227E 

04, 

0.15030E 

01, 

0.13890E 

01, 

0.  13900E 

01. 

0.44800E 

04' 

0.25608E 

04, 

0.15830E 

01. 

0.1*590E 

01, 

0.14590E 

01, 

0.49770E 

04' 

0.27010E 

04, 

0.16620E 

01, 

0.15280E 

01. 

0.1S290E 

01, 

0.54740E 

04' 

0.28434E 

04, 

0.17420E 

01, 

0.15980E 

01, 

0.1S990E 

01, 

0.0 

• 

0.29879E 

0*. 

0.18220E 

01. 

0. 16670E 

01. 

0.16680E 

01, 

0.0 

' 

0.31345E 

04, 

0.19020E 

01, 

0.17370E 

01, 

0.17380E 

01, 

0.0 

' 

0.32834E 

0*. 

0.19820E 

01, 

0.18060E 

01. 

0.18080E 

01, 

0.0 

' 

0.3*343E 

04, 

0.20610E 

01. 

0.18760E 

01. 

0.18780E 

01, 

0.0 

' 

0.35873E 

0*. 

0.21*1 OE 

01. 

0.19460E 

01. 

0.19480E 

01, 

0.0 

' 

0.37424E 

04, 

0.22220E 

01 

DATA  X 2/ 

0.20150E 

01, 

0.20180E 

01, 

0.0 

• 

0.38995E 

04. 

0.23020E 

01, 

0.20720E 

01, 

0.20750E 

01, 

0.0 

t 

0.40586E 

04, 

0.23670E 

01, 

0.20850E 

01, 

0.20880E 

01, 

0.0 

' 

0.42198E 

04, 

0 , 23820E 

01. 

0.215*0E 

01, 

0.21580E 

01, 

0.0 

• 

0.43828E 

04, 

0.2*620E 

01. 

0.22240E 

01. 

0.22290E 

01, 

0.0 

' 

0.45476E 

0*, 

0.25430E 

01. 

0.22940E 

01, 

0.22990E 

01, 

0.0 

' 

0.47144E 

04, 

0.26230E 

Oi, 

0.23630E 

01, 

0.23700E 

01, 

0.0 

• 

0.48829E 

04, 

0.270*0E 

01, 

0.2*330E 

01. 

0,2*41 OE 

01, 

0.0 

* 

0.50532E 

04, 

0.27850E 

01, 

0.25030E 

01, 

0.25120E 

01, 

0.0 

* 

0.52254E 

04. 

0.28660E  01, 

0.25730E 

01, 

0.25830E 

01, 

0.0 

' 

0.53991E 

0*. 

0.29470E 

01, 

0.26A20E 

01, 

0.26550E 

01, 

0.0 

* 

0.55745E 

0*. 

0.30280E 

01, 

0.27120E 

01, 

0.27260E 

01, 

o.o 

* 

0.57515E 

04, 

0.31100E 

01, 

0.27820E 

01, 

0.27980E 

01, 

0,0 

* 

0.59300E 

0*. 

0.31920E 

01, 

0.28520E 

01, 

0.28700E 

01, 

0.0 

• 

0.61 1 0 1 E 

04, 

0.32740E 

01. 

0.29220E 

01, 

0.29430E 

01, 

0.0 

' 

0.62918E 

04, 

0.33560E 

01. 

0.29920E 

01, 

0.30150E 

01, 

0.0 

' 

0.64749E 

04, 

0.34380E 

01, 

0.30620E 

01, 

0.30880E 

01, 

0.0 

• 

0.66607E 

04, 

0.35210E 

01, 

0.31320E 

01, 

0.31610E 

01, 

0.0 

' 

0.68454E 

04, 

0. 360*0E 

01, 

0.32030E 

01, 

0.32340E 

01, 

0.0 

' 

0.70327E 

04, 

0.36870E 

01 

OATA  X 3/ 

• 

0.32730E 

01* 

0.33070E 

01, 

0.0 

' 

0.72215E 

04, 

0.37700E 

01, 

• 

0.33**0E 

01, 

0.3381 OE 

01, 

0.0 

' 

0.7*1 16E 

04, 

0.38540E  01, 

• 

0.3*1*0E 

01, 

0.34550E 

Oil 

0.0 

« 

0.76029E 

0*. 

0.39370E 

01. 

72 
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• 0.34850E  01.  0.35290E  01 

• 0.35560E  01.  0.36040E  01 

• 0.3626 OE  01.  0.36780E  01 < 

• 0.36970E  01.  0.37530E  01 

• 0.37680E  01.  0.38280E  01 

• 0.38400E  01.  0.39040E  01 

• 0.391 10E  01.  0.39800E  01, 

• 0.39830E  01.  0.40550E  01, 

• 0.40540E  01.  0.41320E  01, 

• 0.41260E  01.  0.42080E  01, 

• 0.0  . 0.0 

• 0.70000E-02.  0.68500E-02, 

• 0.70000E-02.  0.70000E-02, 

• 0.70000E-02.  0.70000E-02, 

• 0.69500E-02.  0,70000.^-02, 

• 0.69500E-02.  0.69500E-02, 

0009  OATA  X 4/ 

• 0.69S00E-02.  0.69500E-02, 

• 0.69S00E-02.  0.69S00E-02, 

• 0.69500E-02.  0.69500E-02, 

• 0.69500E-02.  0.70000E-02, 

• 0.70000E-02.  0.70000E-02, 

• 0.70000E-02.  0.70000E-02, 

• 0.69999E-02.  0 .69999E-02, 

• 0.69500E-02,  0.70000E-02, 

• 0.69500E-02.  0.70000E-02, 

• 0.69S00E-02.  0.70000E-02, 

• 0.69500E-02.  0.70000E-02, 

• 0.70000E-02.  0.70000E-02 

• 0.70000E-02.  0.70000E-02, 

• 0.70000E-02.  0.70001E-02 

• 0.69709E-02.  0.6R999E-02, 

• 0.70078E-02.  0.70000E-02 

• 0.69955E-02.  0.70199E-02, 

• 0. 70000E-02.  0.7020 9E-02, 

• 0. 70000E-02.  0.70500E-02 

0010  oat  a x 5/ 

• 0.69999E-02,  0.71000E-02, 

• 0.70000E-02.  0.71000E-02, 

• 0. 70000E-02.  0.71000E-02, 

• 0 , 70000E-02.  0.70999E-02 

• 0.70000E-02.  0.71001E-02, 

• 0.70000E-02.  0.71500E-02, 

• 0 . 70000E-02.  0.71999E-02 

• 0.70000E-02.  0.72000E-02, 

• 0.70000E-02,  0.72001E-02 

• 0. 70000E-02.  0.72500E-02, 

• 0. 70000E-02.  0.73000E-02, 

• 0.70000E-02.  0.73000E-02, 

• 0.7000 IE-02.  0.73000E-02, 

• 0.70500E-02,  0.73000E-02 

• 0.70500E-02.  0.73500E-02 

• 0.70500E-02,  0.74000E-02 

• 0.71000E-02.  0.74000E-02, 

• 0.71000E-02.  0.74000E-02 

• 0 . 71000E-02.  0.74S00E-02, 


DATE  » 75224  13/41/03 


0.0 

, 

0.77954E 

04. 

0.40210E  01. 

0.0 

• 

0.79894E 

04. 

0.41060E  01. 

0.0 

, 

0.81844E 

04, 

0.41900E  01, 

0.0 

• 

0.83807E 

04, 

0.42750E  01. 

0.0 

• 

0.85782E 

04, 

0.43600E  01, 

0.0 

• 

0.87768E 

04, 

0.444S0E  01. 

0.0 

• 

0.89766E 

04, 

0.4S310E  01. 

0.0 

• 

0.91774E 

04, 

0.46170E  01, 

0.0 

• 

0.93796E 

04, 

0.47030E  01. 

0.0 

• 

0.0 

, 

0.47890E  01. 

0.0 

• 

0.0 

, 

0.0  . 

0.48119E 

01, 

0.12S66E 

02, 

0.60000E-02, 

0.49700E 

01, 

0.12606E 

02, 

0.80000E-02, 

0.49700E 

01, 

0.12663E 

02. 

0.80000E-02, 

0.49700E 

01, 

0.12760E 

02, 

0.79500E-02, 

0.49700E 

01, 

0.12865E 

02, 

0.79500E-02  / 

0.49700E 

01, 

0.13003E 

02, 

o.eooooE-02. 

0.49700E 

01, 

0.13163E 

02, 

0.80000E-02, 

0.49700E 

01, 

0.13328E 

02, 

0.79999E-02. 

0.49700E 

01, 

0.13519E 

02, 

0.80000E-02, 

0.49700E 

01, 

0.13716E 

02. 

0.80000E-02, 

0.49700E 

01, 

0.13913E 

02. 

0.80000E-02, 

0.49700E 

01, 

0.14130E 

02. 

0.80000E-02, 

0.49700E 

01, 

0.14343E 

02. 

0.80000E-02, 

0.0 

• 

0.14552E 

02. 

0.80000E-02. 

0.0 

, 

0.14779E 

02. 

0.80000E-02, 

0.0 

• 

0.14987E 

02. 

0.80000E-02. 

0.0 

t 

0.15194E 

02. 

0.79500E-02. 

0.0 

• 

0. 1541 lE 

02. 

0. B0S00E-02, 

0.0 

• 

0.1S610E 

02. 

0.80196E-02, 

0.0 

• 

0.15813E 

02. 

0.79895E-02. 

0.0 

• 

0.16017E 

02. 

0.80000E-C2. 

0.0 

, 

0.16212E 

02. 

0.8046 4E-02. 

0.0 

• 

0.16392E 

02. 

0.80519E-02. 

0.0 

• 

0.16577E 

02. 

0.80500E-02  / 

0.0 

, 

0.16772E 

02. 

0. 81 000E-02. 

0.0 

, 

0,1 694 3E 

02. 

0, 81 OOOE-02, 

0.0 

, 

0.17124E 

02. 

0.61000E-02. 

0.0 

, 

0.1729SE 

02. 

0,81 000E-02. 

0.0 

t 

0.17458E 

02. 

0.81 000E-02, 

0.0 

• 

0.17625E 

02, 

0.81500E-02, 

0.0 

• 

0.17777E 

02. 

0.82000E-02, 

0.0 

• 

0.17932E 

02. 

0.82000E-02. 

0.0 

• 

0.18087E 

02. 

0.82000E-02, 

0.0 

t 

0.1 821 9E 

02. 

0.82001E-02, 

0.0 

• 

0.18469E 

02. 

0, 82SOOE-02, 

0.0 

• 

0.18525E 

02, 

0 . 82999E-02. 

0.0 

t 

0.18581E 

02. 

0.83000E-02. 

0.0 

• 

0.1B830E 

02. 

0.83000E-02, 

0.0 

• 

0.18947E 

02. 

0.83000E-02, 

0.0 

• 

0.19070E 

02. 

0.83S00E-02. 

0.0 

• 

0.19191E 

02, 

0.B3500E-02, 

0.0 

• 

0.19333E 

02. 

0.84S00E-02, 

0.0 

• 

0.19447E 

02. 

0.84500E-02  / 

73 


AEDC-TR-76-15 


f OUTRAN  IV  6 LEVEL  21 


DATE  ■ 7522* 


13/41/03 


0. 

71000E-02. 

0 . 75000E-02, 

0.0  . 

0. 1956 IE  02. 

0. 

71000E-02. 

0.75000E-02. 

0.0  . 

0.19697E  02. 

0. 

71000E-02, 

0.75500E-02, 

0.0  , 

0.1980BE  02. 

0. 

71500E-02. 

0.76000E-02, 

0.0  , 

0.19921E  02. 

0. 

71500E-02. 

0.76000E-02. 

0.0  , 

0.20025E  02. 

0. 

71500E-02, 

0.76000E-02. 

0.0  , 

0.20156E  02. 

0. 

71999E-02. 

0.76000E-02. 

0.0  , 

0.20264E  02. 

0. 

72000E-02. 

0.76000E-02. 

0.0  , 

0.0  t 

0. 

0 , 

0.0  . 

0.0  , 

0.0  . 

0. 

37253E-09. 

0 • 1*90  IE-08. 

0.30928E-05. 

0.30143E-02, 

0. 

Ill 76E-08 , 

0. 1B626E-08. 

0.28610E-07, 

0.3306*E-02. 

-0. 

2*99 8E-0*. 

0.22352E-08, 

0.28610E-07, 

0.73268E-02. 

0. 

15006E-0*. 

-0.24996E-0*. 

0.2B610E-07, 

0.68241E-02. 

-0. 

15009E-0*. 

0.15002E-04.-0.29996E-03. 

0.80565E-02. 

0. 

l*995E-0*. 

-0.15019E-0*. 

0.28610E-07, 

0.99821E-02. 

•0. 

l*991E-0*. 

0.15021E-0*. 

0.19073E-07, 

0.91564E-02. 

0. 

15019E-0*. 

•0. 19988E-0*. 

0.95367E-08, 

0, 10224E-01 . 

-0. 

20021E-0*. 

0.1U76E-08. 

0.2B610E-07, 

0.10937E-01, 

0. 

2*587E-07. 

0,1 1 176E-08.-0.29998E-03. 

0.10333E-01. 

DATA  X 7/ 

0. 

14901E-08. 

-0.29973E-0*. 

0.28610E-07, 

0, 10767E— 0 1 . 

-0. 

24967E-0*. 

0, 1527*E-07, 

0.19073E-07, 

0.1165*E-01. 

0. 

l*992E-0*. 

0.201 17E-07, 

0.0  , 

0. 10010E-01 . 

-0. 

l*991E-0*. 

-0.30012E-0*. 

0.0  . 

0.111 76E-01 . 

0. 

14992E-0*. 

0.2868SE-07, 

0.0  , 

0.10946E-01. 

“0. 

19988E-0*. 

0.0  . 

0.0  , 

0.93161E-02. 

0. 

1 1 176E-08. 

0.0 

0.0  . , 

0.m95E-01, 

0. 

1 1 176E-08. 

0.18999E-0T, 

0.0  , 

0.97370E-02, 

•0. 

27081E-0*. 

-0. 156*62-07, 

0.0  , 

0.96234E-02. 

0. 

39530E-05, 

-0.22055E-05. 

0.0  , 

0.10141E-01. 

0. 

368*2E-0*. 

0.32052E-0*. 

0.0  , 

0.90042E-02. 

•0. 

29083E-0*. 

-0.60603E-05. 

0.0  . 

0.92697E-02, 

0. 

2011TE-07. 

0.20793E-04. 

0.0  , 

0.84900E-02. 

0. 

1*901E-0B. 

-0.19982E-04. 

0.0  , 

0.97B70E-02, 

-0. 

2996*E-0*. 

0.11176E-08, 

0.0  , 

0.B6B6BE-02. 

0. 

1 1 176E-08. 

0.11 176E-0B. 

0.0  » 

0.75607E-02. 

0. 

7*506E-09, 

-0.1825*E-07. 

0.0  . 

0.952*5E-02. 

0, 

26312E-07, 

0.10803E-07. 

0.0  . 

0.66376E-02. 

-0.300122-04.  0 .24994E^>4, 
DATA  X 8/ 

0.28685E-07,-0. 19983E-0*. 
0.0  . 0.15274E-07, 

0.0  . 0.14901E-08, 

0.27940E-07,  0.2A97IE-0*. 
0.0  .-0.2001  IE-0*. 

0.37253E-09.  0.11176E-08. 
0.10431E-07,  0, 7A506E-09. 
0.2*9*3E-0*.  0.37253E-09. 
-0.1A980E-OA.-0.50012E-05. 

0.15005E-0*.  0, 10005E-0*. 
-0.200UE-04*  0.  » ] 1 TOE-00. 
0.11 1 76E- 08.-0. 279*0 E-0 7. 
0.1 1 176E-08.  0 «2*989E-0*. 
•0.300UE-0*.-0,19p  '3E-0*. 
0,11 176E-08.  0 • 1*90  IE-08. 


0. 1956 IE  02.  0.85000E-02. 
0.19697E  02.  0.85000E-02. 
0.19808E  02.  0.85000E-02. 
0.19921E  02.  0.85000E-02. 
0.20025E  02.  0.86000E-02. 
0.20156E  02.  0.66000E-02. 
0.20264E  02.  0.85999E-02, 
0.0  t 0.86001E-02. 

0.0  . 0.0  . 
0.30143E-02,  0.37253E-09. 
0.33064E-02.  0,7*506E-09« 
0.73268E-02.-0.25006E-0*. 
0.682*lE-02.  0.1500AE-0*. 
0. 80565E-02.-0 .200 12E-0*. 
0.99821E-02.  0.25332E-07, 
0.91564E-02.  0.11 176C-08. 
0, 1022AE-01 .-0.29991E-0*. 
0.10937E-01.  0.11176E-08. 
0.10333E-01.  0.11 1 76E-08  / 

0.10767E-01 .-0.29983E-0*. 
0.1165*E-01.  0.111 76E-08. 
0. 10010E-01 . 0. 7*S06E-09. 
0.111 76E-01 . 0.37253E-09. 
0.10946E-01.  0.11176E-0B. 
0.931 61E-02.-0.24985E-0*. 
0.1119SE-01,  0. 49986E-05. 
0.97370E-02.  0.18031E-0*. 
0.96234E-02.-0.28677E-05. 
0.10141E-01.-0.63939E-05. 
0.90042E-02.  0.14984E.03. 
0.92697E-02.-0, l**72E-0*. 
0.84900E-02.  0. 1*61 0E-0*. 
0.97870E-02.-0.20002E-04, 
0.B686BE-02.  0.1U76E-0B. 
0.75607E-02.  0.74506E-09. 
0.952*5E-02.-0.29057E-07. 
0.66376E-02.  0.37253E-09, 
0.8860BE-02.-0.49964C-05  / 


0.70419E-02, 
0.76401E-02. 
0.76996E-02. 
0.964B1E-02. 
0.37BBBE-02.' 
0.27516E-01. 
>0.21B*9E-0  1 .' 
0.20108E-01. 
0.27939E-02, 
0.5A031E-02. 
0.455 94E-02. 
0.53101E-02, 
0.7490SE-02* 
0.35721E-02.' 
0.83603E-02. 


0.99946E-05. 
0.11 176E-0B. 
0.11 176E-08. 
0.20489E-07. 
■0.50206E-05. 
0.10015E-0*. 
-0.85682E-08. 
0.111 76E-08. 
0.U176E-08. 
0,2*985E-0*. 
>0.  15000E-0*. 
0.58079E-05. 
0.14996E-0*. 
■0.200 1 1E-0*. 
0,11 176E-0B. 
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0014 


001S 


0016 


IV  6 LEVEL  21 


HSCP 


DATE  ■ 75224 


13/41/03 


• 0.11176E-08i 

• 0.24986E-04, 
•-0.14997E-04, 

• 0.15015E-04, 
DATA  X 9/ 

•-0.19983E-04 

• 0.10431E-07 
•-0.37253E-08 


-0.50005E-05 

O.1O005E-04 

0.11176E-08 

-0.30006E-04 

0.30028E-04 

-0.17509E-07 

-0.37253E-08 


0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

o.o 


0,465702-02.  0.111762-08, 
0,5580'.E-02.  0.11176E-08. 
0.5766  IE-02.  0.19981E-04, 
0.44830E-02.  0.55879E-08  / 

0.618*4E-02.-0.85682E-08, 
0.252J 3E-02.  0.10431E-07. 
0.0  . 0.14975E-04. 


• 0.0 

•-0.372S3E-10 

•-0.74506E-10 

• 0.14998E-05 
•-0.10004E-05 

• 0 • 1 0006E-05 
— 0.99961E-06 

• 0.99938E-06 
•-0.10012E-05 

• 0.150132-05 
•-0 . 1 490  IE-08 
•-0.55879E-09 

• 0.14982E-05 
•-0.99949E-06 

• 0.99942E-06 
•-0.99949E-06 

DATA  X10/ 

• 0.14991E-05 
•-0.74S06E-10 
•-0.74506E-10 
■ 0 . 1 7086E-05 
•-0.13842E-06 
•-0.14477E.04 

• 0.19540E-0S 
•-0.12666E-08 
•-0.55879E-09 

• 0.19979E-OS 
•-0.74506E-10 
•-0.372S3E-10 
•-0.18999E-08 

• 0.20008E-05 
•-0.18995E-08 

• 0.0 
• 0.0 

•-0.18626E-08 

• 0.0 


0.0 

0.49990E-06 
0.1 1 176E-09 
0.14901E-09 
0.1499BE-05 
0. 1 000  IE-05 
0.10013E-05 
’0.1001 4E-05 
0 .14991 E-05 
0.74506E-10 
0 . 74506E- 1 0 
0. 19979E-05 
’0.55879E-09 
■0.14901E-08 
0.20008E-05 
’0.18999E-0B 


0.0 

0. 15394E-01 
-0.19073E-09 
-0.19073E-09 
-0.19073E-09 
0.19999E-05 
-0.19073E-09 
-0.95367E-10 
-0.95367E-10 
-0. 19073E-09 
0.19999E-05 
•0, 1 9073E-09 
-0. 95367E- 1 0 
0.0 
0.0 
0.0 


0.0  . 0.0  * 

0. 6663 IE -04,-0 .37253E-1 0 » 

0 .29297E -04,-0 .74506E-10. 
0.1 6464E-03.  0.1 5004E-05. 
0.1 060SE-03.-0. 10002E-0S. 
0.75V89E-04.  0.15008E-05. 
0.13184E -03.-0, 15Z74E-08. 
0.60272E-04.-0.5S879E-09, 
0.4715OE-04,  0.19997E-05. 
0.71716E-04.-0.74506E-10. 
0.33874E-04.-0.74506E-10. 
0.77820E-05.  0. 1 9988E-05. 
0.6>902E-04,-0,74S06E-10t 
0.20992E- 04, -0*372532-10. 
0, 11428E-04.-0.37253E-1 0. 
’0,3391 0E-04 .-0 , 74506E-1 0 / 


0.0  , 

0.0 

0.0  . 

0.0 

0.968582-09. 

0.0 

0.67055E-09. 

0.0 

0. 180602-06, 

0.0 

0.960582-05. 

0.0 

0.40721E-06. 

0.0 

0.12892E-05, 

0.0 

0.149882-05, 

0.0 

0.74506E-10, 

0.0 

0.74506E-10, 

0.0 

0.894072-09, 

0.0 

0.74506E-10, 

0.0 

0. 14999E-05. 

0.0 

0,149842-05. 

0.0 

0.55879E-09, 

0.0 

0.37253E-09, 

0.0 

0,1 4986E-05, 

0.0 

0. 15007E-05, 

0.0 

• 0,7141 lE-04.  0.1 4991 E-05. 
.-0.23804E-04.-0.74506E-10. 

. 0 . 14343E-04.-0. 13033E-05. 

, 0.35095E-04,  0.90897E-07, 

. 0,3051 8E-05.  0.57567E-06. 

. 0,491 33E-04.-0.49470E-04, 
.-0.17853E-04.  0.98299E-06. 

. 0.52643E-04.-0.98027E-06, 

. -0. J35692-05.  0. 150032-05. 
.-0.85449E-05.-0.74506E-10. 

. 0. ^9335E-04,-0.37253E-l  0. 
.-0.67291E-04.  0, 19372E-08. 

. 0.10040E-03.-0.37253E-10. 
.-0.31891E-04,  0.49967E-06, 

♦ 0.3601 IE-04.. 0.499752- 06. 

♦ 0 .6561 3E-0 5.-0 , 7450 6E- 10. 

• 0.36621 E-05.-0.74506E- 10. 
.-0.20294E-03.-0 . 1 0431E-08. 

. 0.58166E-03.  0.50113E-06  / 


DATA  XU/ 
•-0.372532-10 
•-0.74506E-10 
•-0.14967E-05 

• 0.99860E-06 
•-0. 1 00032-05 

• 0, 15007E-05 
»»0. 745062- 10 
•-0.74506E-10 

• 0.20007E-05 
•-0. 745062-10 
•-0.74506E-10 


-0*745062-1 0< 
-0.37253E-10’ 
-0.37253E-10’ 
0.5001 6E-06’ 
-0.S0027E-06’ 
-0.74506E-10 
0. 18626E-08 
•0.  14994E-05’ 
0. 1 4988E-05 
•0.1 1 176E-09 
0.50012E-06 


O.C 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 

0.0 


..0. 164932-02. 
0. 16461E-02. 
-0.51270E-03. 
0.20386E-03. 
0.50049E-04. 
-0.34790E-04. 
0*-21002. 
-0.1 1978E-03. 

0.141142-03. 
•0,1 0452E-03. 
0.19662E-04. 


0*501242-06* 
0.89407E-09. 
0.74506E-10. 
0.74606E-10. 
0.14990E-05. 
0.100012-05, 
.558792-09*  . 
0. 999792-06, 
0.150072-05* 
0,74506e-10, 
0.745062-10, 
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-0.14991E-05 
0.99979E-06 
-0.10010E-05 
0 « 14984E*05 
*0. 74506E*  1 0 
0.40978E-09 
0.0 

OATA  X12/ 
0.38170E  02 
0.42992E  02 
0.45770E  02 
0.45813E  02 
0.47B18E  02 
0.49386E  02 
0.50685E  02 
0.51806E  02 
0.52798E  02 
0.53692E  02 
0.0 
0.0 
0.0 
0.0 

0.5S178E-01 
0.3234OE-O1 
0.23941E-01 
0.21529E-0 1 
0.17679E-01 
OATA  *13/ 
0.13769E-01 
0.11787E-01 
0.10378E-01 
0.93121E-02 
0.85026E-02 
0.0 
0.0 
0.0 
0.0 

-0.40361E-04 
-0.37794E-04 
0.17000E-03 
-0.58597E-05 
-0.20860E-04 
-0.35399E-05 
-0.32197E-05 
-0.30831E-05 
-0.30690E-05 
-0.31018E-05 
DATA  X14/ 

0.0 

0.0 

0.0 

0.0 

-0.59220E-06 

-0.33915E-07 

-0.29616E-03 

-0.89275E-07 


>0 .50027  £*06 
■0.74506E-10 
0.20003E-05 
•0 .200 1C|E*05 
0.1341  ,E*08 
0.4097 JE-09 
0.0 


0.413S5E  02 
0.462'. 8E  02 
0.49014E  02 
0.49047E  02 
0.51091E  02 
0.52722E  0? 
0.54198E  02 
0.55  ?97E  02 
0.56161E  02 
0.57320E  02 
0.0 
0.0 
0.0 
0.0 

0.S8152E-01 

0.32470E-01 

O.23878E-01 

0.21688E-01 

O.18278E-01 


0.14525E-01 

0.12603E-01 

0.11143E-01 

0.10001E-01 

0.91142E-02 

0.0 

0.0 

0.0 

0.0 

-0.40848E-04 
-0.37315E-04 
0.15203E-03 
>0, 33331 E»05 
-0.21507E-04 
-0.37267E-05 
-0.37805E-05 
-0.366B5E*05 
-0.34713E-05 
-0.35350E-05 


0.0 

0.0 

0.0 

0.0 

■0.58375E-06 

■0.43852E-07 

■0.26809E-03 

■0.91442E-07 


0.71400E  01 
0.71800E  01 
0.90500E  01 
0.10640E  02 
0.12080E  02 
0.13280E  02 
0.14280E  02 
0.15170E  02 
0.15980E  02 
0.16720E  02 
0.17410E  02 
0.18070E  02 
0.18700E  02 
0.0 

0. 2308 1E-0 1 
0.20129E-01 
0. 16B50E-01 
0.15092E-01 
0.13371E-01 


0.10629E-01 

0.92143E-02 

0.84111E-O2 

0.76691E-02 

0.70499E-02 

0.67125E-02 

0.64501E-02 

0.61500E-02 

0.0 

■0.77265E-03 
>0.  10091E-04 
-0.10925E-04 
■0.35603E-05 
-0.13714E-04 
■0.47140E-05 
-0.13971E-05 
-0.19137E-05 
>0,1 8855E-05 
>0.1  1237E-05 


■0.75072E-06 

■0.15002E-05 

■0.15013E-05 

0.0 

-0 . 89R76E-05 
■0.42040E-07 
0.14243E-07 
■0.33627E-07 


0.54932E-05 

•0.36316E-04 

0.13763E-03 

-O.52940E.O4 

0.86670E-04 

0.0 

0.0 


0.42758E  02 
0.47769E  02 
0.51072E  02 
0.51127E  02 
0.53830E  02 
0.56122E  02 
0.58126E  02 
0. 5991 OE  02 
0.61522E  02 
0.62992E  02 
0.0 
0.0 
0.0 
0.0 

0.58339E-01 

0.36818E-O1 

0.3C295F-01 

0.28413E.01 

0.25040E-01 


0.21044E-01 

0.18663E-O1 

0.16795E-01 

0.15278E-01 

0.14056E-01 

0.0 

0.0 

0.0 

0.0 

-0.31659E-04 
•0,30206E*04 
0.1154OE-03 
>0 . 77686E*05 
>0 . 23648E-04 
>0.63094E*05 
«0,59986E*05 
■0.50714E-05 
■0.51110E-05 
■0.481 02E-05 


0.0 

0.0 

0.0 

0.0 

■0.50628E-06 

■0.20882C-07 

•0.22392E-03 

•0.60635E-07 


■0.74506E* 

■0.99860E* 

■0.52154E* 

0.40978E- 

•0.74506E* 

■0.49889E* 

0,0 


0.38396E 

0.41916E 

0.45106E 

0.45155E 

O.47404E 

0.49334E 

0.50891E 

0.52249E 

0.53464E 

0.54570E 

0.0 

0.0 

0.0 

0.0 

0.66549E- 

0.37253E- 

0.27234E- 

0.24710E- 

0.207R9E* 


0.16430E-01, 
0.14233E-01, 
0.1 ?656E*0 1 • 
0.11463E-01, 
0.10586E-01, 
0.0  , 

0.0  , 

0.0  . 

0.0  , 

■0.47519E-04, 
-0.43164E-04. 
0.15985E-03, 
■0.33868E-05. 
•0.25080E.04, 
■0. 38276E»05t 
■0.38131E-05, 
•0.32553E-05, 
•0.33221E-05. 
■0,311 18E-05  / 


0.0  , 

0.0  t 

0.0  . 

0.0  . 

■0.65973E-06. 
■0.53501E-07, 
>0. 30344E-03* 
-0.10811E-06. 
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13/41/03 


0020 


0021 


0.R7354E-0B. 

0 . 1 8?88£»07 . 

0.0  • 

0.24433E-07. 

0.21892E-07. 

-0.42454E-07, 

•0. 39231 E-07.-0.1S716E-07, 

•0.37324E-07, 

•0.4772SE-07, 

-0.255P6E-07. 

•0.23469E-07, 

-0.17457E-07, 

-0.22272E-07. 

-0.2713AE-07, 

•0.14973E-07, 

•0.13582E-07. 

•0.11 977E»07 » 

•0.16753E-07. 

-0.18055E-07, 

•0.65234E-08, 

-0.64328E-08. 

-0.80682E-08, 

•0.66679E.08, 

-0.71093E-0B, 

-0.42401E-08, 

•0.40766E-08, 

-0.37S51E-08, 

-0.54486E-08, 

-0.54412E.08, 

0.0  • 

0.0  . 

-0.37443E-08. 

0.0  , 

0.0  , 

0.0  • 

0.0  . 

-0.74506E-12. 

0.0  , 

0.0  , 

0.0  • 

0.0  . 

0.29802E-11. 

0.0  , 

0.0  , 

0.0  • 
OAT*  *15/ 

0.0  . 

0.0  • 

0,0  , 

0.0  / 

0.69580E  01. 

0.6981 0E  01. 

0. 7961 OE  01. 

0.696 1 OE  01. 

0.77340E  01. 

0.79690E  01. 

0.70200E  01. 

0.88740E  01. 

O.802S0E  01, 

0.70230E  01, 

0.88960E  01. 

0.B0270E  01. 

0.71960E  01, 

0.98770E  01. 

0.81860E  01, 

0.74310E  01. 

0.10666E  02. 

0.84150E  01. 

0.767O0E  01, 

O.U310E  02. 

O.86760E  01. 

0.78630E  01. 

0.11846E  02, 

0.89S40E  01, 

0.80630E  01. 

0.12293E  02. 

0.92460E  01. 

0.B2120E  01, 

0.12667E  02, 

0.9SA70E  01. 

0.0  • 

0.0  • 

0.0  . 

-0.25556E-03. 

0.54B7SE-02, 

-0.16S27E-03. 

0.23493E-03, 

0.11356E-01, 

0.32042E-03. 

0.14S90E-02. 

O.11798E-01. 

0.82102E-03, 

0.16891E-02, 

0.11629E-01. 

0.129S6E-02, 

0.21 R28E-02, 

0.86783E-02, 

0.2019RE-02, 

0.23782E-02, 

0.69620E-02. 

0.2SA75E-02. 

0.23619E-02, 

0.57708E-02, 

0.27283E-02, 

0.19795E-02. 

0.48289E-02, 

0.28715E-02. 

0 . 1 6236E-02 , 

0.40369E-02, 

0.29941E-02, 

0.13516E-02. 

0.34077E-02. 

0 • 30325E-02 » 

0.0  , 

0.0  • 

0.0  . 

0.36619E-05, 

0.25903E-05. 

0.25014E-05. 

-0.127BAE-05, 

C.34HSE-05, 

0.25451 E"05. 

0.13946E-03, 

0.24256E-03, 

0.16496E-03. 

-0.37O97E-0S. 

-0.25071E-04. 

0.1 5B99E-05, 

0.30618E-05. 

-0.64868E-05, 

0.28290E-05. 

0.51733E-06, 

-0.37476E-05, 

0.67945E-07, 

-0.31326E-05, 

-0.29OA8E-05. 

0.11945E-06, 

-0.1B257E-05, 

-0.2B476E-05. 

0.22802E-06. 

-0.12885E-0S, 
0ATA  *16/ 

-0.261A1E-0S, 

0.92646E-07, 

-0.98218E-06, 

•0.25038E-05  / 

-0.3492SE-10, 

0.0  . 

0,0  . 

0.0  , 

-0.B0632E-0B. 

0.17834E-06. 

-0.48627E-09, 

0.51037E-07, 

-0.78710E-08. 

0,3461 0E-10 , 

-0.2784BE-04, 

-0.10383E-03. 

-0.1322AE.06, 

0.41190E-07, 

0.687T1E*07, 

0.13541E-07. 

-0.13900E-07, 

-0.13966E-07, 

-0.127UE-08. 

-0.39916E-08, 

-0.1A723E-07, 

0.SS730E-0B, 

0.B1368E-08, 

-0.12029E-07, 

0.39T94E-08, 

0.30976E-09, 

-0.7M82E.0B, 

0.2S660E-08, 

-0.47721E-09, 

•0.3S454E-08, 

0.66161E-09, 

-0.13392E-OB, 

•0.27310E-08, 

0.74990E-09, 

0.0  , 

0.0  . 

0.0  / 

0022 

T»TX/1 .8 

0023 

GO  TO  (10. 

0024 

1000 

H.-l 

0025 

GO  TO  50 

0026 

10 

GO  TO  (11, 

0027 

20 

GO  TO  (13, 

0028 

1003 

H«-2 

0029 

GO  TO  50 

0030 

30 

GO  TO  (31 

0031 

31 

H*  -2 

0032 

GO  TO  50 

0033 

11 

IF  ( IT  ,LT 

0034 

T 1 ■ 100. 

0035 

T2»  110. 0 

0036 

INC*  10 

0037 

IF  (K  ,EQ, 

0038 

NK*  52 

(31  • 13<31 « 13, 1 3) ,K 


(T  ,GT,  600.1)  GO  TO  IS 
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0039 

GO  TO  25 

0040 

12  IF  l(T  .LT.  298. 

IS)  .OR. 

0041 

TI«  298.15 

0042 

T2a  300.0 

0043 

iNCa  100 

0044 

NKa  14 

0045 

GO  TO  25 

0046 

13  IF  (IT  ,LT.  100. 

) .OR.  IT 

0047 

TIa  100.0 

0048 

T2a  200.0 

0049 

INCa  100 

0050 

NKa  11 

0051 

25  NK2a  NK-2 

0052 

00  28  Ila  1.NK2 

0053 

T3«  T ? • I NC 

0054 

IF  (T3  .EO.  300. 

1 T 3*  298, 

0055 

IF  ( (T3  .EO.  308, 

>15)  .OR. 

0056 

IF  (IT  ,GE.  Tl)  , 

>ANO.  IT  , 

0057 

TIa  T2 

0058 

T2a  T3 

0059 

28  CONTINUE 

0060 

Ila  NK-1 

0061 

GO  TO  40 

0062 

26  NKa  51 

0063 

NK2a  NK-2 

0064 

00  27  IIal,NK2 

0065 

T3a  T2« INC 

0066 

IF  (IT  ,GE.  Tl)  . 

ANO.  IT  . 

0067 

TIa  T2 

0068 

T 2a  T3 

0069 

27  CONTINUE 

0070 

0071 

0072 

0073 

0074 

0075 

0076 

0077 

0078 
0070 
0080 
0081 
0082 

0083 

0084 

0085 

0086 

0087 

0088 
0080 
0000 
0001 


DATE  • 75224 


• GT.  1500.1)  GO  TO  15 


.GT.  1000. 1 1 GO  TO  15 


•LT.  T21)  60  TO  40 


40 
1001 

41 


42 

43 


IS 

50 


Ila  NK-1 

GO  TO  141,42,43),  KK 
H«-l 

GO  TO  50 
TC«  T-Tl 

if  ik^eq.'”  Return  CI"<,n,3,,#TC  * ciik.ii,2))*tc  • ciik 

H»  , 1 8E04*M 
RETURN 
TCa  T-TI 

RETURN2<K,II,4I#TC  * C^,,(,,,,3,  ,#TC  * C2IK.II.2)  )*TC  ♦ C2IK 

TCa  T-Tl 
Ka  K-2 

IF  IK  ,EQ.  0)  Ka] 

U*  ' IC3IK,II,4)*TC  ♦ C3 IK , 1 1 , 3) ) *TC  ♦ C3IK, II ,2) ) *TC  ♦ C3IK 

IF  IK  .EG.  3)  Ka; 

RETURN 
Ha  o.o 
RETURN 
END 
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FORTRAN 

0001 

0003 

0003 

0004 


0005 


0006 


IV  0 LEVEL  31 


COEFF 


OATE  « 75334  13/41/03 


C 


r 


c 


SUBROUTINE  COEFF (COOE .ARRAY ) 

OIMENSION  ARRAY  ( 731 ,HHA(73> .HI (73) .51 (73) .CPI (73) *HH1 I73i . 
1 HM3(73) ,HW3(7Z) .HHCI40) 

INTEGER  COOE 

CHEMICAL  • SENSIBLE  ENTHALPHY  FOR  T » TI. 

OATA  mhI/ 


19.3937814E.4 

, 4.9679389 

. 0.0 

« 

0.0 

• H 

31 .0650436E  «S 

. 4.9819393 

. 0.0 

• 

0.0 

, 0 

31.7343991E*4 

. 6.5117735. 

3.4740577E-4 

t 

0.0 

, OH 

43. 781 0593E  »4 

, 7.7703187, 

1.3914575E-4 

t 

0.0 

. NO 

5*4. 995035 7E  *4 

, 7.4991057. 

1.6384351E-4 

« 

0.0 

, CO 

64.3795667E*? 

. 6.3305090, 

3. 7333736E-4 

• 

0.0 

• H3 

7-330.5134204, 

6.8101303.  3. 

8505033E-4  , 0. 

0 

, 

N3 

8 -630.6330969 

. 7.3851083 

. 3.5993393E-4 

, 0.0 

. 03 

9-. 0615 

, 4.9683153 

. 0.0 

8 

0.0 

, AR 

A-175 045. 83 47470.  13.506037 

1 0.0 

• 

0.0 

, C03 

8-. 0615 

. 4.9683153 

. 0.0 

• 

0.0 

. NE 

C-103633. 3476190,  7.1088046 

. 7.678847E-4 

% 

0.0 

• H30 

0-1304.1576077 

. 7.6133733 

. 1 . 36571 70E-4 

,0.0 

• N3 

E-1340. 9978438 

• 7.8839888 

, 1.5887864E-4 

,0.0 

• 03 

F -0.03187 

, 4,9681371 

• o.o 

• 

0.0 

• AR 

0-177175.4038470.  14.394766 

. 0.0 

• 

0.0 

. C03 

H -0.03187 

• 4.9681371 

1 0.0 

• 

0.0 

. NE 

1-104643.9857996.  8.5658373 

• 5. 041 8784E-4 

.0.0 

/ H30 

ENTROPY  FOR  T 

> TI. 

DATA  HH3/ 


1 

36.367556 

. 5. 377978 IE-03 ,-0.9431 089E-06, 

0.7157348E-10, 

H 

3 

37.535609 

, 5.420 3863E-03 .-0 • 9540  069E-06 , 

0.7257386E-10, 

0 

3 

48.450387 

. 5. 341 4646E-03.-0.41 77444E-06 . 

0.0  , 

OH 

4 

48.070143 

. 0.8479758E-03.-1.4094686E-06. 

1.0545863E-10, 

NO 

5 

45.057500 

. 0.8187987E-03.-1.3406314E-'*, 

u»9s.'5*B43E-lO* 

CO 

6 

39.869398 

. 0.7339605E-03.-1 . 30 1 9850E-06 « 

0.9363449E-10, 

H3 

7 

57,375034 

, -1.8300467E-03,  1 . 7559B3OE-05 

. -0.581 5B33E-08. 

N3 

8 

47.335630 

. 0. 7577403E-03.-0. 8381 940E-06, 

0.0  . 

02 

9 

36.446448 

. 4, 55710] 4E-0 3. -0.511 9033E-06. 

0.0  , 

AR 

A 

47,518398 

. 1 .1513786E-03.-1 .3066643E-06, 

0.0  • 

C02 

8 

34.4116313 

. 4.5570565E-03.  -0.5118949E-06 

■ 0.0  , 

NE 

C 

43.361707 

. 0. 834488 3E -03.-0 , 7763385E-06 

,0.0  . 

H20 

0 

46.856474 

. 5. 0506993E-03.-0. 3700634E-06. 

0.0  • 

N3 

E 

50.3049115 

. 5. 3780447E-03.-0. 3853719E-06, 

0.0 

02 

F 

38.5134376 

. 3.0363745E-03.-0.3383395BE-6. 

0.0  , 

AR 

G 

51.5035705 

• B.44473E-03  .-0.61701 16E-06, 

0.0  , 

C02 

H 

34.8197333 

• 4.5503533E-03.-0.6903733E-06, 

0.0463530E-09, 

NC 

I 

45.6963147 

. 0.6353374E-03,-0.*115987E-06. 

0,0  / 

H20 

SPCCFIC  HEAT 

FOR  T > TI, 

OAT»  Hr,;/ 

l 

4,98615 

,0.0  ,0.0  , 

0.0  , 

H 

3 

5,0648570 

,-4.98170  •"  1 E— 5 . 0.7133534E-08. 

0,0  , 

0 

3 

6.1708494 

. 0.60407t"E-3.  0.4777381E-0T, 

•1  .448839 3E -11. 

OH 

4 

5.909456 

• 1 .9345802E-3,  -0, 45604 16E-06, 

0.3931985E-10. 

NO 

5 

5.5733319 

• 3.00766511-3,  -0.4553B45E-06, 

0.3TT70T6E-10, 

CO 

6 

9,6700579 

.-4, 367650 5E -3.  3.53645S4E-06. 

•0.5547857E-09, 

H2 

7 

5.7766563 

. 1.5061552E-3.  -3.0759936E-07. 

0.0  , 

N3 

8 

5,607584 

• 3.7614501E-3,  -0.88881 1 3E-06. 

1.0954859E-10. 

02 

9 

4,96815 

.0,0  , 0,0  , 

0.0  , 

A* 

A 

6.680396 

. 0.5990060E-3,  -1 . 7096379E-06. 

0.1809910E-09, 

C02 

79 
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OATE  « 7522* 


13/*l/03 


0007 


0008 


0009 


8 *.96815 
C 7,1 864387 
6.5727317 
7.6813079 
*.96*15 
10.8521173 
*.96815 


0.0  . 0.0  . 0.0 

0.9855705E-3.  0 ,*6*0690E-06, - 1 .0002320E-1 0 . 
0.9o99083E-3.-0.9596*9*E-07  , 0.0 
*.*2812*SE-*,-l.R738208£-08 
0.0  . 0.0 

1.6107530E-3.  -1 . 7 1 25528E-07 , 0.0 

0.0  . 0.0  . 0.0 

2.5861179E-3.-2.3692669E-07,  0.0  / 

SENSIBLE  ENTHALPY  COEFFICIENTS  FOR  T .L£.  T 1 .900R-I 8009 , 
18009-2 700H. 2700R-*000». 


0.0 

0.0 


1 5.976*766 


OATA  HI/ 
•272.57872 
•-*9.9253121 
• 0.19368 


• 0.19368 


•-.0615 


.0615 


• -0.02187 


-0.02187 

-1858. *935816.  8.5657282.  S.0*20885E-* 
ENTROPY  COEFFICIENTS  FOB  T 
2700H-*0009  REVISED  EQ. 

OATA  51/ 

0.0155587  , 

1.5770590E-02. 

• 1 .291 8283E-P?  , 

2. 0*7 1 31 6E-02, 

. 1 .291 70*2E-02. 

1.8092773E-02,- 
-1.8300*67E-02, 
0.7577*02£-02, 

. *.557101*E-03. 

1.1512786E-02, 

. *.S570S65E-03. 
0.8?**883E-02. 
5.0506993E-P3. 
5.2T80**7£-03. 


0.0 


.LE.T 1,9009-1 8009. 18009-27009. 


• 39.25139* 

• *2.269726 
•31.5883510 

• *1.78*58 

• 29.5537807 

• 37. *37166 

• 57,27502* 

• *7.225630 
•36. **6**80 

• *7,518398 
•3*. *116312 

• *3.261707 

• *6.856*7*0 

• 50.20*9115 
•38.512*276 

• 51,5035705 
•3*. 8197332 

• *5.69631*7 


R,***72E-03  , - 

*.550352?E-03. 
0.635327AE-O2,  . 
SPECIFIC  HEAT  COEFFICIENTS  FOR 
18C09-27009. 27009**0009  REV, 

Data  cpi/ 

0 .-1.301 5396E-03 , 1 

1 , 2.2*l5270E-O3,-C 

0,0  , 0. 


• 7,23667810 

• 5.78*6285 

• *.96815 


80 


NE 

H20 

N2 

02 

AR 

C02 

NE 

H?0 


. 6.28*9888  . *.2*60662E-*  , 0.0 

• 

N2 

. 6.7005911  , *.60979*1E-*  , 0.0 

• 

02 

. *.9682152  , 0.0 

. 0.0 

• 

AR 

6.5377190,  28. 320 1 S 7£  — * 

. -0 , 3858033E -6 

t 

C02 

• *.9682152  , 0.0 

, 0.0 

• 

NE 

.913576.  8, 1 3*9063E-*  . 

0.0  . 

H20 

6,8101302.  2.8505033E-* 

• 0.0  , 

N2 

7.3851082.  2 . 5992293E-* 

,0.0  , 

02 

. *.9681*18  , 0.0 

• 0.0 

t 

AR 

11.078**,  5.*166925E-* 

. 0.0 

• 

C02 

• *.9681*18  , 0.0 

. 0.0 

• 

NE 

.1089372,  7.6785528E-*  . 

o 

o 

• 

H20 

7.6123722,  1.3657170E-* 

, 0.0 

• 

N2 

7,8839888,  1 .588786*£-* 

. 0.0 

• 

02 

, *.9681271  . 0.0 

, 0,0 

• 

AR 

12.7*7913.  2.31 06225E-* 

, 0.0 

• 

C02 

. *.9681271  , 0.0 

, 0.0 

• 

NE 

H20 


0.5615997E-05 

0.93376*E-09 

* 

N2 

0.539O*27E-05 

0.8573795E-09 

• 

02 

-5 « 3921 089E-06 

0,966391 7E-09 

• 

AR 

0. 63*  71 25E-05 

0.9526361E-09 

* 

C02 

-5,3907*1 0E-06 

0.9659557F-09 

• 

NE 

.6362909E-05 

1 ,07*831 7E-09 

* 

H20 

. 7559820E-05 

-0.5815R32F-08, 

N2 

0.8281 9*E-06  , 

0.0 

* 

02 

-0.51 19033E-06 

0.0 

• 

AR 

1 .2fl666*3E-06 

0,0 

* 

C02 

-0,511 89*9£-06 

0.0 

» 

N£ 

0.7782285E-06 

0.0 

• 

N20 

0 . 370063*E-06 

0.0 

• 

N2 

3. 385271 9E-06 

0.0 

• 

02 

-0,2283? 95 8E-06.  0.0 

• 

AR 

0.6170) 16E-06 

0.0 

« 

C02 

-0,690273 3£ -06, 

0,0*6?520E-09 

• 

NE 

0.*1 15987E-06  , 0.0 
T.LE.T1 ,9009-18009, 

/ 

H20 

.5759008E-06,  -0. 37*6857£-09 

• 

N2 

,*58fl023E-06» 

.0 

• 

02 

.0  . 

.0 

• 
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FORTRAN  lv 


0010 


0011 

0012 

0013 

001*. 

0015 

001ft 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 

0038 


G LEVEL  21 


COEFF  DATE  * 75224  13/41/03 


C 


5.3131429 

4.96815 

7.5482768 

5.7766562 

5.60768400 

4.96815 

6.6803960 

4.96815 

7.1564387 

6.5727317 

7.6813079 

4.96815 

10.8521173 

4.96815 

5.9764766 


0.8363827E" 

0.0 

4.2815697E" 
1.50615S2E- 
2.761450 1E« 
0.0 

0.5990060E- 

0.0 

0.9855705E* 

0.9099083E* 

4.4281245E- 

0.0 

1 , 61 07530E« 

0.0 

2.5861 179E‘ 


02.-0. 
. 0. 
04.  0. 
03.-2. 
03,-0. 
, 0 

02,-1. 
. 0 

03,  0. 

03, -0. 

04, -1. 

. 0. 

03.-1. 
. 0. 
03.-2. 


3088999E-05 
0 

7157717E-06 
0759936E-07 
88881 13E-06 
0 

7096379E-06 

0 

464069E-06. 

9596494E-07 

8738208E-08 

0 

7125528E-07 

0 

3692669E-07 


• 0,44931 04E-09 

• 0.0 

• -1.349S489E-10 

. 0.0 

, 1.0954859E-10 

• 0,0 

. 0, 1 8099 1 OE-09 

• 0.0 

-1 .0002320E-1 0 
• 0.0 
• 0,0 
0,0 
0.0 
0.0 
0.0 


G1885  FREE  ENERGY. 


DATA  MHC/ 

1 5. 3150524E*04, -15. 6253583 

2 6. 0670223E*04, -21. 8518477 

3 1 • 08438 12E *04, -24. 82 13338 

4 2. 2896290E»04. -28. 1298227 
S-2.5226419E*04,-26.3634922 

6 ] ,44l5362E*03. -17. 8897838 

7 1.2087000£*03. -25. 6082461 

8 1,281 489 3E-0  3.-27 .3792772 
9-9. 2608408E»04. -27. 957881 1 
A-S.6574736E •04,-25.0072600 


.-9.9767035E-03.  0.5610564E-67  , 
.-0.9782325E-03,  0.56141 14E-07  , 
.-1.6149044E-03,  0, 7304  7 18E-0 7 • 
.-l  • 597936 1E-0 3,0, 8680 263E-0 7 , 
.-1.5649177E-03.  0.8457005E-07  , 
.-1.3549595E-03,  0.6781 319E-07  , 
.-1.5325905E-03.  0.8179650E-07  . 
•-1.6271943E-03,  0.8783605E-07  , 
.-2.5501816E-03.  1 .3506145E-07  , 
.-1.8841889E-03.  0.8839907E-07/ 


00  1 1*1 ,72 

1 ARRAY  < 1 ) * 0.0 

GO  70(10,20,30.40,60,60,60,80,90) .CODE 

10  DO  11  1*1.40 

11  ARRAY ( 1 1 * HHC ( 1 ) 

RETURN 

20  00  21  1*1,72 

21.  ARRAY  ( 1 ) « HM1  (1) 

RETURN 

30  00  31  1*1.72 

31  ARRAY ( 1 1 ■ MM2 11) 

RETURN 

40  DO  41  1*1,72 

41  ARRAY ( 1 ) ■ MM  3(1) 

RETURN 

50  00  51  1*1,72 

51  ARRAY ( 1 ) * MMA ( l ) 

RETURN 

60  00  61  1*1,72 

61  ARRAY(l)  ■ Ml (It 
RETURN 

80  00  81  1*1,72 

81  ARRAY ( 1 ) ■ 51(1) 

RETURN 

90  DO  91  1*1,72 

91  ARRAY ( 1 ) * CPI (1) 
return 


END 


C02 

NE 

R20 

N2 

02 

AR 

C02 

NE 

H2D 

N2 

02 

AR 

C02 

NE 

M20 


M 

0 

OH 

NO 

CO 

H2 

N2 

02 

C02 

M20 
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FORTRAN 


0001 


000? 

0003 

0004 

0005 

0006 
0007 

ooos 

0009 

0010 
0011 
0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 
0022 

0023 

0024 

0025 

0026 

0027 

0028 

0029 

0030 

0031 

0032 

0033 

0034 

0035 

0036 

0037 


FORTRAN 

0001 


0002 

0003 

0004 

0005 

0006 
0007 


IV  G LEVEL  21  CCC  OATE  » 75224  13/41/03 

function  cccim.n.t) 

C TEMPERATURE  T IS  IN  OE GHEES  R 

C N2  -1,10.15,  02-2,20,25,  AR-3,30,35,  C02-4,40,45,  NE-5,50,55 
GO  TO(100.200.3C0) , M 
100  GO  T0(1,2,3,4,5),N 

C 3H0  VIR1AL  COEFFICIENT (Cl - AEDC-  TR-71-39. 

1 CCC  » <3.5528210  -4.7253632E-02  •T)/(1.0  -4.7446815E-03  »T> 

GO  TO  1000 

? CCC*  ( -4. 234S758E-01 I / ( 1 , 0 -0.9728086E-02  »T) 

GO  TO  10C0 

3 CCC  »(((-l. 0787834E-07  »T)  ♦ 2.4336141E-04)  • T - 

1 0.19016416)  • T • 59.879538 

GO  TO  1000 

4 CCC  * 115.17385  -0.1338398  • T • 0.3190R38E-4  • T • T 
GO  TO  1000 

5 CCC  *0.0 
GO  TO  1000 

200  GO  TO(10,20,30,40,50),N 
C T (OC/OT ) 

10  CCC  ■(-30.396628E-03  «T) / < < (22.5120025E-06  »T) -9.4893630E-03) 

1 * T . 1.0) 

GO  TO  1000 

20  CCC»(4.1194317E-03  «T)/(1.0  - 0.9728086E-02  »T) 

GO  TO  1000 

30  CCC* ( < ( (-3.2363502E-07  • T)  • 4.8672282E-04)  • T 
1 -0.19016416)  • T) 

GO  TO  1000 

40  CCC  * 0.1338398  * T ♦ 0.6381676E-4  • T • T 
GO  TO  1000 
50  CCC  * 0.0 
GO  TO  1000 

300  GO  TO(15.25.35.45.5S) ,N 
C T«»(D»*C/OT**) 

15  CCC  * -28.8444637E-05  • T • T ♦ 136.B577931E-0B  • T »»3/ 

1 ( 1.0-9. 4893630E-3  * T ♦ 22.5120025E-6  * T »T)  ••  2 
GO  TO  1000 

25  CCC  * (-4.0074185E-05  • T * T)  / (1.0  -0.972B086E-0?  »T) 

GO  TO  1000 

35  CCC  * 4.8672282E-04  »T  »T  -6.4727004E-07  »T  »T  »T 
GO  TO  1000 

45  CCC  * 0.6381676E-4  • T • T 
GO  TO  1000 
55  CCC  ■ 0.0 

iooo  ccc  ■ ccc  ♦ ioo.o 

RETURN 

END 


IV  G LEVEL  21  SONV  OATE  ■ 75224  13/41/03 

SUBROUTINE  SONV(VTM.F .T.CP.G.A) 

C TEMPERATURE  T MUST  BE  IN  DEGREES  R 
C SUBROUTINE  TO  CALCULATE  SONIC  VELOCITY 
IF (F.EO.O) VTM*2B«9646 
G ■ CP/ (CP- ( 1 ,987165/RTM) ) 

X ■ 32.16  • 1.987165  • 1054,3502645/  1.355B1T9 

A * (ABS((G»X»T)/VTMI) 4*0.5 

RETURN 

ENO 
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FORTRAN 


0001 


0002 

0003 

0004 

0005 

0006 

0007 

0008 


0009 

0010 

0011 

0012 

0013 

0014 

0015 

0016 

0017 

0018 

0019 

0020 
0021 


0022 

0023 

0024 

0025 

0026 

0027 

0028 


0029 

0030 

0031 

0032 


0033 

0034 

0035 

0036 


IV  G LEVEL  21  8BB  DATE  • 75224  13/41/03 

FUNCTION  B88 (M.N, T ) 

C TEMPERATURE  T IS  IN  OEGREES  R 

C N2  -1,10,15.  02-2,20,25,  AR-3,30,35,  CO2-4.40.45,  NE-5,50,55 
GO  TO( 100 ,200 ,300 ) ,M 
100  GO  TO C 1 . 2, 3. 4, 51 ,N 

C 2NO  VIRUL  COEFFICIENT  (B)  - AEOC  - TR-71-39. 

1 BBB  * (((-0.3534783E-03  • T)  - 0.6636132)  * T • S.  1364637E*02) 

1 / (((-0.9296798E-05)  • T - 1.7158299E-02)  * T ♦ 1.0) 

GO  TO  1000 

2 BBB  * ( ( (-1.5050447E-04  *T>  - 4.7667818E-01 ) * T ♦ 4.2BS7o84E*02> 

1 / ( ( (-0 . 5095962E-05*  T) • 1 .41230 1 0E-02)  * T ♦ 1.0) 

GO  TO  1000 

3 8BB  • ((((-1.4860291E-08  • T)  ♦0.47T2262E-03)  • T - 

1 0.8705351)  • T • 3. 8904656E  *02)  / < (( <1 .4973286E-10  • 

2 T)  ♦ 1.1790085E-05)  • T -1 .4845725E-02)  * T ♦ 1.0) 

GO  TO  1000 

4 88R  * (((  0.7026469E-04  • T)  -5.0336R43E-01 ) * T *5.31S2906E*02) 

1 /(((-2.2018674E-0B  *T>  -0 .61 1 35B4E-02) *T  *1.0) 

GO  TO  1000 

5 BBB  *(59.425695-0.27144987  *T) / 1 1 . 0-1 .6463468E-02  *T> 

GO  TO  1000 

200  GO  T0(10.20.30.40,50) ,N 
C T (OB/DT) 

10  BBS  *(«3.4560T92£*02)  / ( ( <-0.29627746E-05*T)  -1 . 1286423E-02) 

1 • T « 1.0) 

GO  TO  1000 

20  BBB  *(-3.1691728E*02>/(  ( ( t 0.4514126E-09  *T)  -0 . 3228332E-05) 

1 *T  -0,9521 19BE-02)  • T *1.0) 

GO  TO  1000 

30  888  «(-2.T524532£*02  » 4.4273617E-02  • T)  / ( 1 .0  - 
1 0.94B551 IE-02  *T) 

GO  TO  1000 

40  BBB  *(((((  18.9162693E-12  »T>  -7.878B948E-09)  • T 

1 -0.9807759E-03)  • T • 5.3458176)  • T) / ( ( ( (-7. 1340504E-08 

2 *T ) -1 , 104452E-02)  • T • 1,0)**  2) 

GO  TO  1000 

50  888*  (-6, 19928  TOE*  01  • 1 .8B86040E-02  *T) / ( 1 .0-1 . 7014989E-02  *T) 

GO  TO  1000 

300  GO  T0(15. 25. 35*45.55) »N 
C T**(D**B/OT**) 

15  BBB  ■<6.7503388E*02>/( ( ( (1.3S34608E-09  *T)  -0.4815325E-05) 

1 *T  -0.8951382E-02)  • T ♦ 1.0) 

GO  TO  1000 

25  BBS  •(((((  0.3089713E-09  »T)  -O.’-.BBSOSTE-OS)  *T  *l,093566TE-02) 

1 *T  -4,9627303) *T  -3.3796998  ♦01)/(((« (0. 9466770E-10 

2 *TI  -0.2105136E-06)  *T  ♦ 1 .2893683E-04)  • T 

3 -1.9686162E-02)  * T ♦ 1.0) 

GO  TO  1000 

3S  880  ■ ( 3.9858659E*02  -3.4251675E-02*T) /( 1 .0*0.661 1836E-02  *T) 

GO  TO  1000 

45  BBB  -(<:(< (-50.1 159356E-16*T) -47. 4547 667E-13) *T 

1 - 1 ,446291 7E -07) * T-12S.9326564E-05)  * T 

2 *11,8074171) *T*T) / ( ( ( (-7, 1340S040E-06  *T) 

3 -1.104452E-02)  • T ♦ 1.0)  **  4) 

GO  TO  1000 

55  8B8* 1 1 ,58520796*02-3, 1 3T8493E-02  *T)/(1.0-1.8918656E-02  *T) 

1000  RETURN 
Em n 
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h- 

i' 


NOMENCLATURE 


A* 
a -i. 

B* 

B 

C 

CP 

CP 

f 

H 

h 

MW 

m 

n 

P 

R 

S 

s 

T 

V 

XM 

XN 

n 

SUBSCRIPTS 


(See  Eq.  (19)) 

Coefficients  for  equations 
(See  Eq.  (20)) 

Second  virial  coefficient,  cm® /mole 

Third  virial  coefficient,  cm® /mole 

Specific  heat  at  constant  pressure,  Btu/lbm-mole-°R 

Specific  heat  at  constant  pressure,  Btu/lbm-°R 

Fuel-to-air  ratio,  lbm^ue^/lbma^r 

Enthalpy,  Btu/lbm-mole 

Enthalpy,  Btu/lbm 

Molecular  weight,  lbm/lbm-mole 

Fuel,  hydrogen-to-carbon  weight  ratio,  lbmjj/lbmQ 

Mole  fraction,  moles/ mole 

Pressure,  psia 

Gas  constant 

Entropy,  Btu/lbm-mole°R 
Entropy,  Btu/lbm°R 
Temperature,  °R 
Specific  volume,  cm® 

Number  of  hydrogen  atoms  in  a mole  of  hydrogen  fuel 

Number  of  carbon  atoms  in  a mole  of  hydrocarbon 
fuel 

Efficiency 


a Pressure  in  atmospheres 

B Base  temperature 

eq  Equilibrium 
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I 

i 

j,k,m,n 

P 

T 

u 


Value  where  the  effects  of  chemical  dissociation  are 
first  considered 

Constituent 

Property  coefficients 

Pressure 

Temperature,  °R 

Universal 
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